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^! Abstract: We construct a theory of fields living on continuous geometries with frac- 

^ ! tional Hausdorff and spectral dimensions, focussing on a fiat background analogous 

to Minkowski spacetime. After reviewing the properties of fractional spaces with 
fixed dimension, presented in a companion paper, we generalize to a multi-fractional 
scenario inspired by multi-fractal geometry, where the dimension changes with the 



Q ! scale. This is related to the renormalization group properties of fractional field the- 

*^ \ ories, illustrated by the example of a scalar field. Depending on the symmetries of 

O 



the Lagrangian, one can define two models. In one of them, the effective dimension 
flows from 2 in the ultraviolet (UV) and geometry constrains the infrared limit to 
be four-dimensional. At the UV critical value, the model is rendered power-counting 
renormalizable. However, this is not the most fundamental regime. Compelling ar- 
rS \ guments of fractal geometry require an extension of the fractional action measure to 

5^ \ complex order. In doing so, we obtain a hierarchy of scales characterizing different 

geometric regimes. At very small scales, discrete symmetries emerge and the notion 
of a continuous spacetime begins to blur, until one reaches a fundamental scale and 
an ultra-microscopic fractal structure. This fine hierarchy of geometries has implica- 
tions for non-commutative theories and discrete quantum gravity. In the latter case, 
the present model can be viewed as a top-down realization of a quantum-discrete to 
classical-continuum transition. 
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1. Introduction 

With the birth of quantum gravity models, it became clear that ordinary geometry 
is inadequate to describe the microscopic texture of spacetime. Using the concept 
of spectral dimension, borrowed from the spectral theory of fractal geometry, it was 
realized that many scenarios are characterized by a scale-dependent dimension. The 
spectral dimension ds is an indicator of the effective number of directions a pointwise 
probe feels when diffusing in a given ambient spacetime for a short amount of time. 
For ordinary manifolds, it corresponds to the integer topological dimension D, but 
for fractals this may not be the case. Several quantum gravity or quantum space- 
time scenarios such as causal dynamical triangulations (CDT) ||l], §, asymptotically 
safe quantum Einstein gravity [^, spin foams 0-0, Hofava-Lifshitz gravity 0, §], 
K-Minkowski non- commutative fleld theory |^, and non-local super-renormalizable 
quantum gravity |[10| are deflned in ambient spacetimes with D = 4 dimensions, but 
the spectral dimension at small scales differs from that value and ds < D (in many 
cases, ds ~ 2) in the UV. The change of dimensionality with the scale is another 
typical property of fractals (more precisely, multi-fractals). Thus, it may be natural 
to regard all these models as different manifestations of the fact that the application 
of quantum mechanics to spacetime itself leads, in general, to a fractal geometry. 

It is important to establish whether this is only an analogy or not, because the 
reduction of dimensionality is intimately related to the renormalization properties 
(meant as ultraviolet flniteness) of quantum gravity. In particular, there seems to 
be a conspiracy between UV flniteness and a spectral dimension (is ~ 2 at small 
scales |Tl|-|[13|. However, quantum gravity research has not fully exploited the vast 



fleld of fractal geometry, and the spectral dimension alone is insufficient both to 
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characterize a physical process as "fractal" and to control its geometric properties. 
First, other notions of dimension can and, actually, must be compared with rfg in 
order to classify the geometry more precisely. Second, finding the spectral dimension 
at a given scale is a far cry from having dimensional flow under full control at all 
scales. Third, there exist other details in the geometric and topological structure 
of fractals which were never or seldom checked for in quantum gravity. In order to 
better understand the connection between quantum gravitational physics and fractal 
geometry, one can take two opposite perspectives. In one, a given model of quantum 
gravity can be chosen and its fractal properties checked. Unfortunately, due to 
technical difficulties, in the great majority of the cases the spectral dimension is the 
only computable fractal indicator. In the other perspective, one could start from 
fractal geometry itself and attempt to construct a theory of quantum gravity with 
the desired properties (dimensional flow from 2 to 4, UV finiteness, and so on). 

The second approach was advocated and outlined in its main qualitative features 
r5|. To have a much closer contact with fractal geometry, a more rigorous 



m 

programme has been initiated in [^, |T7|. Exploiting the characterization of many 
fractals as systems governed by a fractional differential structure, we constructed 
a continuous space with Euclidean signature and whose Hausdorff and spectral di- 
mensions are non-integer. This empty space, called fractional Euclidean space and 
denoted by £^ , has a notion of distance, volume and dimension. In particular, it 
has topological dimension D and Hausdorff and spectral dimensions d-^ = d^ = Da 
(for non-anomalous diffusion), where < a < 1 is a fixed real parameter. Further- 
more, it is endowed with symmetries, although not the usual rotation and translation 
group. In practice, on fractional Euclidean space one can ask the same questions and 
perform all the operations allowed in ordinary Euclidean space, but with a different 
calculus. 

In this paper, we build upon the results of ll6| and carry out the agenda spelled 



out therein in greater detail. The idea is (i) to extend £^ to a space with Lorentzian 
signature, (ii) realize dimensional flow via tools of multi-fractal geometry, (iii) discuss 
an example of fleld theory and its renormalization properties, and (iv) generalize to 
models with even more realistic fractal properties. It will turn out that step (ii), 
which was not really made in [|T3|-[|T3|, contains some pleasant news: The renor- 



malization group (RG) flow and the multi-fractal construction are one and the same 
entity, but described with two different languages. This may be unsurprising, since 
the RG flow is based upon a scale hierarchy just like multi-fractals. What is per- 
haps surprising is the quantitative match: once symmetries are given, purely fractal 
geometric considerations lead to the same total action prescribed by an almost- 
traditional fleld theory analysis. On top of that, geometric requirements can flx the 
effective dimension of spacetime in the infrared (IR) to c/h = <is ~ 4, provided the 
dimension in the UV is 2. Also, step (iv) will be crucial to probe scales even smaller 
than those at which the RG flow takes place, and will allow us to make contact with 
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some features discovered numerically in causal dynamical triangulations. The over- 
all physical picture has been shortly presented in [|17|. Gravity is not included, yet, 
but fractional Minkowski spacetime will be sufficient to illustrate the basic features 
of the proposal. We leave the fractional extension of general relativity to another 
publication. Below is the plan of the paper. 



Section ^. Following the Euclidean construction of \TE\, we define a fractional 
ambient space with Lorentzian signature and describe its geometry (section 
ZTD and its Hausdorff and spectral dimensions (section |2.2|) . The choice of 



coordinate presentation and the role of the boundary, which were not discussed 
in ||16i, are here treated in detail. The fractional generalization of Lorentz 



transformations is presented in section 2.3 



Section |^. Spacetimes with a multi-fractal structure are introduced. Exploiting 
the lore of multi-fractal geometry, the measure in the action is argued to be a 
linear superposition of contributions with fixed dimensionality (sections |3.1| and 
Ol) . The dimension of spacetime changes with the scale, but not arbitrarily: 



the dimensions in the UV and in the IR are, in fact, intimately related (section 
p.3|) . In particular, we shall analyze the role of two and four dimensions at, 
respectively, small and large scales. Observational constraints on the fiow of 



the dimension near the infrared limit are discussed in section 3.4. 



Section ^. An example of classical fractional field theory is provided by a real 
scalar field 0, whose ultraviolet finiteness can be easily probed via a standard 
power-counting argument (section [4.11) . Generic classical fractional actions are 
constructed in section [4.2| , where the equation of motion of is also derived. 
The choice of integer or fractional Lorentz symmetries in the Lagrangian density 
lead to two independent models, respectively, the integer symmetry scenario 
(section [4.3|) and the fractional symmetry scenario (section [4.4| ). The Green 
function inverting the kinetic operator is calculated in section [4.5| . The power- 
counting renormalizability of the models is discussed at several points in the 
section, in particular in section W^. 



Section |^. To get in closer touch with fractal geometry, fractional theories are 
extended to the case where the measure is a linear superposition of fractional 
measures of complex order. Combining this superposition in a real quantity, one 
obtains a measure with logarithmic oscillations (section |5.1| ). The average of 
the measure over a log-period corresponds to the real-order fractional measure. 
The oscillations are due to a discrete symmetry (sections |5.2| and ^.3|) , and 
require a redefinition of the Hausdorff and spectral dimension in line with the 
definitions employed in fractal geometry (section |5.4|) . 
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Section [^. After summarizing the physical picture in section |6.1| , we outhne 
a research agenda focussed on the quantum theory, the inclusion of gravity 
and cosmological applications (section |0| ). Section ^]3| is devoted to connec- 
tions with doubly special relativity, non- commutative spacetimes and quan- 
tum gravity approaches. Multi-fractional field theories can be regarded either 
as stand-alone models of quantum gravity or as effective descriptions of other 
theories in certain regimes. In the bulk of the paper we assume the first atti- 
tude. Considering instead the second case, we advance possible applications of 
multi-fractional geometries to other, independent models of quantum gravity 
in section |6.3.3| . 



1.1 Comparison with early proposals 

All the present material is novel with possibly two minor exceptions. The first is 



section |^, which makes heavy use of the results of [0. These, however, are here 
immediately extended to an ambient spacetime with Lorentzian signature. The sec- 
ond exception is the very idea that spacetime be fractal. Before embarking ourselves 
in the construction of multi-fractional spaces with dimensional flow, it is useful to 
draw an exhaustive comparison between our proposal and other spacetime models 
in non-integer dimension which appeared in the early literature. This comparison, 
which was premature for the abstract Euclidean fractional space of fixed dimension- 



aly presented in |I6[ , can give the reader a bird's eye view of the state of the art of 
fractal spacetimes and the status of our theory within. 

Most of previous "fractal" field theoretical proposals assumed spacetime to have 
a non-integer but constant, non-dynamical dimensionality, thus fixing the attention 
to 4 — e dimensions with < e ^ 1 (low-lacunarity regime). Here we mainly 
focus on the rather few papers attempting to realize the texture of spacetime via a 
mathematical fractal or fractal-inspired construction. 

• Quantum field theories in non-integer dimension were considered in [|1^ via a 
simple analytically continued integration. They were regarded as abstract mod- 
els created for the purpose of shedding some light into four-dimensional field 
theory, and the non-integer dimension formally appearing in manipulations was 
not associated with a fundamental modification of spacetime geometry. The 
hope was to unravel properties of the renormalization group flow which do 
not depend on dimensionality. At variance with these dimensionally-continued 
field theories, the philosophy of later publications [^-[^ was that spacetime 
has a non-integer but fixed, non-dynamical dimensionality. The attention was 
focussed on (in = 4 — e for obvious empirical reasons. In other words, these 
models share many aspects with dimensional regularization but the parameter 
e is taken to be physical and non-vanishing, albeit small, (i) In ||I9| a mathe- 
matical justification to dimensional continuation was given, with applications 
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to statistical mechanics and field theory. An axiomatic description of a metric 
space with non-integer dimension was proposed. Neither is this a vector space 
nor is it embedded in a vector space of integer dimension D, so it can hardly 
be compared with our framework. The lack of a manageable set of natural 



coordinates did not allow the author of |T^ to explore the model much, and 
objects such as a "Pythagorean" coordinate distance, volumes, and the symme- 
try group of the space were not investigated. Also, this construction becomes 
more and more complicated with the increase of the topological dimension [^ . 
Yet, it was possible to define an invertible Fourier transform (an open problem 
in fractal geometry, but solved in fractional spaces [0) and a natural Laplacian 



operator, (ii) In pO[, particle physics was defined directly on sets with general 



Borel measure of fixed dimensionality, i.e., a Lebesgue-Stieltjes measure with 
possibly fractal support. Renormalization properties in electrodynamics were 
considered for a general measure with low lacunarity, dn = A — e. Convergence 
of the Feynman diagrams is better than in four dimensions, as it was checked by 
looking at their superficial degree of divergence. This was presented as a "new 
regularization method," meaning that physical applications should be sought 



out only at low lacunarity. (iii) In pT| , |22[| , scalar- field theory in Euclidean 



signature was constructed on low-lacunarity Sierpinski carpets with Hausdorff 
dimension dn = 4 — e. Because of the explicit fractal construction, the scaling 
property of the measure is discrete and the system displays a certain symmetry, 
called discrete scale invariance, which we shall discuss in section |5.2| . The prop- 
agator on a fractal lattice was computed using techniques which would have 
been later developed for the spectral theory on fractals |^. (iv) The interest- 
ing UV properties of field theories in fractional spacetimes and the breaking of 
parity and time reversal therein, all topics we shall amply discuss, were also 



appreciated in |p5[|, where it was noticed that fractional and curvature effects 
are similar in the limit of almost-integer dimension. Even in this case, attention 
was limited to (in = 4 — e.^ 

In contrast with all these approaches, field theories on a genuinely and deeply 
"anomalous" spacetime {dn much different from 4) have received less attention, 
despite their promising applications in modern cosmology and quantum gravity. 
In [^, it was observed that the fine-scale structure of a quantum mechanical 
particle path is very irregular and described by a nowhere-differentiable curve. ^ 
Ord drew inspiration from this fact to propose a model of quantum mechanics in 



^Breaking of time reversal and non-conservation of probability |14 in fractional quantum me- 
chanics were discussed also for fractional generalizations of the Schrodinger equation pS, |2^-~p0{. 

^The Hausdorff dimension of the path is 2 in a classical spacetime g^ and smaller than 2 
in quantum spacetimes with a minimal length ||33| . In the latter case, dfj can be even negative, 
corresponding to a Planckian regime where quantum fluctuations of the spacetime texture are large 
and the particle path is an empty set. 
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fractal spacetime [^]. The idea that spacetime be a mult i- dimensional fractal 



was also hinted at in |]35[. Physical descriptions and implications of fractal 



spacetimes have begun to be focussed on a geometrical perspective in p3|-p!5 
within a Lebesgue-Stieltjes approach. Away from the low-lacunarity regime, 
which misses all the potentialities of fractal spacetimes, the Lebesgue-Stieltjes 
formalism [^ |TB[ can say little unless one specializes to specific measures. 
General absolutely continuous measures away from c^h ~ 4 were discussed in 



[]14| , but their special status as naive "fractal" measures did not allow us to 



make much progress with an adequate level of rigorousness. 



In |T3| , |n|, fractional calculus served as a motivation for the introduction of 
general Lebesgue-Stieltjes models of spacetime, but it was not used in its full 
power. The framework of these papers shares only the main qualitative char- 
acteristic with the present one (namely, anomalous scaling of the measure and 
better renormalization properties of field theories) and it does not possess the 
richness of physical implications we can appreciate here. In particular, while in 
the models of |jl3|]-|jl5[ it is not clear how to construct a rigorous definition of 



an invertible unitary transform between configuration and momentum space. 



this is can be done for fractional spaces [^ thanks to the factorization of the 
coordinates in the measure. 

Through the notion of distance and the calculations of ball volumes, we saw in 
[|16| and will see in section |2| that geometric coordinates q can be interpreted as a 



coordinate system intrinsic to the fractional space, and the mapping x -^ q{x) 
relates the embedding (or extrinsic) viewpoint to a geometric (or intrinsic) 
viewpoint. These two equivalent pictures were conjectured, without giving de- 
tails, at the end of [Q . In |TB|, |I1| , a study of the deformed Poincare algebra of a 
Lebesgue-Stieltjes model with arbitrary measure reached the same conclusion, 
naming the embedding and geometric pictures, respectively, conservative (as a 
system made of two dissipative parts) and dissipative (as a system dissipating 
energy-momentum in a bulk). 

When this article was being finalized, the author became aware of an inde- 



pendent proposal for a fractal spacetime |^-[^. While ordinary spacetime 
is described by smooth differentiable manifolds, abandoning the requirement 
of differentiability naturally leads to a fractal geometry characterized by non- 
absolute scales which can only be measured relatively to one another. Starting 
point, motivation, terminology and tools greatly differ from ours. Central is the 
so-called principle of scale relativity. Motivated by fractal geometry, coordinate 
frames are made explicitly dependent on the scale and, crucially, scales must 
transform according to certain very natural laws generalizing ordinary contrac- 
tions/dilations. This principle can be applied to a number of systems, not only 
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in physics, leading to a severe modification of our perception of Nature. In tlie 
realm of physics, no formal theory has been constructed upon scale relativity, 
but some of the consequences of this principle presents intriguing similarities 
with our fractional approach, including a breaking of parity symmetry and 
an almost obvious dimensional flow. A careful comparison between fractional 
spacetime formalism and the scale relativity proposal may be of mutual ben- 
efit. For instance, in the external scale picture of section ^TT] we gave a few 
a{i) profiles as toy examples realizing a running from some critical value a* 
to the integer charge a = 1. The external scale i is the very same continuous 
scale labelling Nottale's fractal coordinates, and by the simple scaling argu- 
ments of scale relativity one can prescribe the scale dependence of the fractal 
dimension. In our notation, these fractal coordinates are indeed the geometric 



variables q'^{£) = q^[a{£)] and eq. (14) of [^ reads 



When i ^ i^:, a{i) — )■ 1, while at small scales a{i) — )■ a*. On one hand, 
implementation of scale relativity arguments could fix some loose points of 
fractional spacetime models and sharpen the overall physical interpretation. 
On the other hand, fractional models could provide the missing theoretical 
framework wherein to embed the scale relativity principle.^ 

The dependence on the scale (or resolution) can be implemented at the level 
of fields rather than coordinates. Such is the philosophy of wavelet field theory 



41, ^ which, not surprisingly, has better UV properties than ordinary field 
theory. We believe that also this approach converges to the same physics of 
fractional and fractal models. 

Dimensional flow may be realized also in non-fractal scenarios, (i) Field the- 
ory actions with exotic measures were the subject of |^. There, spacetime 



was described by a set of "continuous" coordinates and its dimension formally 
constrained by a variational principle. Despite the fractal-inspired motivation, 
there is no obvious point of contact with fractal geometric scenarios. The 



cosmological "decrumpling" model of [|J] also considered dimension to be a 



dynamical field, (ii) Another unrelated appearance of an exotic Hausdorff di- 
mension in field theory is in four-dimensional gravity with a quantum conformal 



factor [^. As a trace anomaly effect, the conformally invariant IR fixed point 



of this model is associated with an anomalous dimension greater than the topo- 



logical one. (iii) The universe described in [|U|-[|S| has a crystal-like "layered" 



^Such a theoretical framework was also proposed in |Q, where non-differentiable manifolds 
were defined in a rather abstract fashion. Physical applications and consequences of these "fractal 
manifolds," and the connection with standard fractal-geometry tools, are presently unclear to us. 



structure governed by a hierarchy of scales along the topological directions. Due 
to the fact that transitions from one dimension to another are rather sharp, 
away from the transition points these models strongly resemble dimensionally 



regularized spacetimes (see also [^). As we have just seen, and contrary to 
what advertized in [^, ^, the paradigm that the effective dimensionality of 
spacetime depends on the probed scale is far from being new. At any rate, 
fractional spacetime theory is quite different from the crystal-world proposal 
except at the IR fixed point. In particular, in our framework there is no geo- 
metric reason why events in a (in = 2 regime should be planar (two-dimensional 
fractals may not be embedded in a plane: Brownian motion is an example), 
and gravitational waves should be produced even in (in = 3 dimensions, mainly 
for the reasons advanced in [pO|].^ 

The measured value of the dimension of spacetime may be slightly smaller 
than four due to quantum fluctuations and to the intrinsic finite resolution of 
experiments. Then, the infinitesimal covering in the definition of Hausdorff di- 
mension cannot be physically realized in the real world. This effect was studied 
in 1^, where no assumption was made on the true dimensionality of space- 
time, and it is unrelated from field theories living in fractal geometry. Heuristic 
finite-resolution effects of quantum fluctuations have been also considered in 
The change of dimension in the multi-fractal flow considered here is far 



more dramatic than this type of corrections, which can be safely ignored. 

To the best of our knowledge, the following features have been developed here for the 
first time: (i) the construction of a fractal-like structure in Lorentzian signature; (ii) 
an explicit realization of dimensional flow in explicitly multi-fractal spacetime struc- 
tures; (iii) the construction of a system whose UV and IR dimensionality are deeply 
related to each other rather than being fixed phenomenologically; (iv) a detailed clar- 
ification of what we mean by "fractal spacetime" and the fractal interpretation of 
fractional spacetime models; (v) a description of the symmetries underlying fractional 
spacetimes under the perspective of fractal geometry, and their consequent identifi- 
cation with the isometry group in the language of field theory; (vi) the construction 
and physical interpretation of fractional spacetimes with oscillatory measures and 



the associated emergence of discrete scales. Point (iii) was vaguely foreseen in [|14 



(i) and (v) were initiated in [|T3|, |^; the first part of (iv) was completed in |T6[. Point 
(vi) also opens up a novel connection with non-commutative spacetimes, as we shall 
see below. 



Comparing with the proposal in |T3|-|T5[, there are many differences and an 
appreciable amount of novelty. There, the simplified setting of an absolutely con- 
tinuous measure captured the main qualitative features of the central idea, but its 

^To this list of non-fractal references, we should also add the effective field-theory fractional 
equations of motion of |5lt|-[p3. These models seem to be purely mathematical. 
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level of rigor was not satisfactory. On the other hand, here we use the formalism of 
fractional calculus, through which we control all the elements of the proposal more 
strictly and in greater detail. The results of |T3|, [l^ partially lie on a Weyl-type 



integral, a special case of fractional integral. Whenever a point of contact is possible, 
we shall compare the present model with those results and clarify how issues of the 
latter are solved here. 

We stress that, while the great majority of the works cited above considered 
dimensional flow only as a general concept or paradigm, our aim here and in |]16| is 
to construct a hands-on theory of spacetime and fields where one can have as much 
control as possible. The concreteness of the model and its close contact with the lore 
of fractal geometry will allow us to go well beyond the isolated feature of anomalous 
scaling. 

2. Fractional Minkowski spacetime 

The extension of fractal geometry to an ambient spacetime with Lorentzian signature 



is, to the best of our knowledge, a topic virtually untouched in the literature [14 



Nevertheless, after having defined fractional geometry in Euclidean space |T6[, we 
can move to a Lorentzian spacetime straightforwardly. Curvature is not considered. 

2.1 Definition 

We define fractional Minkowski spacetime A4^ oi order a as a D-dimensional em- 
bedding Minkowski spacetime M^ endowed with a set of rules Calc" = {d°', /", d} 
of integro-differential calculus (symbols denote derivative, integration and external 
differential), a measure ga with a given support, a natural norm || ■ ||, and a Laplacian 
/C: 




(2.1) 



2.1.1 Embedding and calculus 

The embedding M^ is Minkowski spacetime in D topological dimensions, with 
"mostly plus" signature (—,+,■■■ , -|-). The embedding coordinates x^ are labeled by 
Greek indices /z, z/, . . . running from to D — 1. The time direction will be sometimes 
denoted as t = x°. 

The embedding is, actually, only a choice of metric once a differential and metric 
structures are defined. The rules Calc" of differential geometry are given by fractional 
calculus [57, Q. This is reviewed in [16, section 2]; here we recall the main definitions 



in one dimension. Given a real coordinate variable x defined on an interval [a;o,a;i], 
the left fractional integral of order a of a function f{x) is 

1 /"'' r/r' 

(^"/)W^=fr^/ , ,1-J i^')^ (2-2) 

r(a) Jxo [x - x'y " 
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where 

0<a<l. (2.3) 

Similarly, the right fractional integral of order a is 

1 r^ dr' 

{rf){x) ' 



r(a) j^ [x' - xy ° 



where integration now is from x to the end of the interval. Also, the left and right 
Caputo derivatives of order a are defined as 

{dV){x) := (/^-"9/)(x) 

1 /"^ dr' 

= . / , , d,,f{x') , (2.5) 

r(l-a)4 {x-xT ^ ^ ^ ^ 

{d"f){x) := {p--df){x) 

1 r^ dr' 

= -TTT^ — 7 / ^-^^5,'/(a;') , (2.6) 

r(l-a) A (x'-x)- ""■"• ^' ^ ^ 

where d is the ordinary first-order partial derivative. Sometimes we will indicate the 
integration terminals explicitly as subscripts, d" = d^^ ^ and 9" = d"^._^. Under the 
transformation 

X ^ Xq + Xi — X , (2.7) 

the left operators are mapped into right operators: 

(r/)(x) = {rF){xo + x,-x), {d^f){x) = (9"F)(xo + x^ - x) , (2.8) 

where F{x) := /(xq + Xi — x). 

One can extend a to other ranges and define other differential operators (e.g., the 
Riemann-Liouville derivative). Different sets of fractional operators can correspond 



to inequivalent fractional spacetimes. In [16| we have justified the use of the Caputo 



derivative via several arguments; the main one is that differential geometry and 
tensor calculus are considerably simplified. 

When xq = — oo, the left operators are called Liouville differintegrals, while for 
Xi = +CXD they are called Weyl differintegrals; they will be all denoted by a subscript 
oo to the left of the symbol. When regarded as an approximation in the limit t ^ to, 
the Liouville operator is employed in mechanics to describe "steady state" systems, 
that is, systems which evolved well after the initial transient phase at to- To get the 
Weyl differintegral from the Liouville differintegral, it is sufficient to set Xq = —Xi in 
(p.8|) and then take the limit Xq — )■ — oo: 



Urf)ix) = WF)i~x) , iooB'^mx) = {ood"F){-x) , (2.9) 

where F{x) := f{—x). 
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Some examples of fractional derivatives and integrals are the following. The 
Caputo, Liouville and Weyl derivatives of a constant are zero, 

a°l = a°l = oo5"l = oo^"l = , (2.10) 

while the left derivatives of a power law are 

d''{x-x,f = -^^^^^ix-x,f-\ /3^0, (2.11) 

i (p — « + Ij 

„s"(x - ^.f = (-1)-° ''!r\ ■ fy\i '^ - -•>'"" • p.12) 

T[p + 1 — a) sm[7r(p — a)\ 

for any x*. The second expression is ill defined for /3 = a. Otherwise, it is real under 
certain conditions on the values of a and /3 and the sign of x — x*; consistently, it 
vanished for /3 = 0. The a-th order left integrals of power laws are given by the 
analytic continuation of the above formulae for a — )■ —a for any /3: 

r(^-^.f = r(''^:;\) (^-^.r°. (2.13) 

„r(i-i.)'' = tlir£(zp^(^_^.)s+<.. (2.14) 

The eigenfunctions of the Liouville derivative operator are exponentials, 

^aV^ = A"e^% (2.15) 



but those of Caputo derivatives with finite Xq are Mittag-Leffler functions |59 



d'^E^[\{x - xqY] = XE4X{x - xo)"] . (2.16) 

Analogous formula can be obtained for the right derivative by making use of eqs. (|2 
and ( p.9| ). For example, the right version of eqs. ( |2.11| ), ( p.l2| ) and ( p.l5| ) are 

a"(a:,-x)^= ^.^^^ + ^\. (xi-x)^-", (3^0, (2.17) 

T{(3 - a + 1) 

oqU [X — Xif) = _ , „ - r— : ; r-z —[X — X^fj , [Z.io) 

T[p + 1 — a) sm[7r(p — a)\ 



9V^ = (-A)V^ (2.19) 



2.1.2 Measure 



The measure and the integration range of the action can be fixed once and for all by 
some simple arguments, which we develop in one dimension. 

The action will be defined via the right integral /". To extend the integration 
range as much as possible while keeping the fractional measure well defined, the 
extremum x must remain finite while taking the limit Xi — t- +oo. The fractional 
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space has then a boundary at some finite x. Without any loss of generahty, we can 
set X = in the action integraL In fact, a translation x -^ y = x + x^^ changes the 
coordinate presentation of fractional operators, but it does not change the physics. 
For instance, the one- dimensional Weyl integral of a function becomes 



J"/ 



1 

1 



dx' 



{x' 



X 



il— a 



fix') 



+ 00 



dx' 



oo 



^ {x' + X* 

dy' 



y 



ll—a 



yy-a 

fiy' + y-x 



fix') 



In the new reference frame, the coordinate dependence of / and of the boundary 
have changed, but this only modifies the way the fractional one-dimensional space 
is embedded in R. As an embedding coordinate change, this operation is always 
possible and one can further set y = x*. For the Liouville integral the measure range 
is in the negative semi-axis, but this would be yet another change of presentation; 
so we can just pick the Weyl integral. The combination oo-^" + oo-^" would hide the 
presence of the singularity at x = by a fictitious integration on the whole line. 

In D topological dimensions, each embedding direction x^ is associated with a 
"fractional charge" a^ (the subscript n is not a vector index). The simplest case is 



of an isotropic fractional spacetime, where all a^ are equal [|T3|, |TJ, |T^. In particular, 
the time direction t = x*^ is on an equal footing with spatial coordinates, cto = a. 
In another configuration worth mentioning for its applications in non-commutative 
spacetimes ([|r^, |BD[ and section |6.3.2| ), time is an integer direction, cto = 1, in which 
case the spatial part of the measure carries the whole effect of dimensional fiow. 
Keeping ao general, the action in fractional Minkowski space is 



(2.20) 





f+oo 


r+oo 


r+oo 


Sa = 


/ A^ = 


/ dga{x)C = 


/ d^XVaix) C, 




Jo 


Jo 


Jo 



where £ is a Lagrangian density, the integral is 

Jo „_n JO 



fi=0 
1 



r(ao; JO 
and the measure along each direction is 



D-l 



+°° dt 

fl-ao 11 Y(a) ./o X 



M=l 



dx 

l-a ' 



(2.21) 
(2.22) 



q" ■■= Qaixn ■-- 



X 



tj.\a 



T(a + 1] 



x^> 0. 



(2.23) 
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To keep notation light, we wrote a instead of a^ in eqs. (|2.20|) and (|2.23|) , with the 
understanding that the time coordinate may scale differently. 

While the measure of Minkowski spacetime is the ordinary Lebesgue measure 
q{x) = ®/_jX^, A^^ is equipped with the Lebesgue-Stieltjes measure 



D-l 

Qaix) := Qaoit) (^ Qaix'') , 

whose scaling property is 

^„(Ax) = A"°+(^-^)"^(x), A>0, 



(2.24) 



(2.25) 



In eq. ( p.23D we have written the measure in the /i direction as a coordinate 
q'^, dubbed geometric or fractional in []TB|. Depending on the differential calculus 
associated with the action integral, {q'^} can be regarded as the natural coordinate 
system spanning fractional spacetime [|l^, section 3.3]. For instance, the fractional 
fundamental theorem of calculus states that the Caputo derivative is the left inverse 
of the fractional integral defined on the same sector (left or right) and the same 
interval ||16|, section 2.3.3]. If we based the calculus on A^^ only on one sector, the 
right integral would be associated with the Weyl derivative 









go 



QaqfM 



T(l-a) 



dx'^ 



xi^ 



ix'^ - x^'Y 



d,. 



(no sum over yu) 



(2.26a) 
Because of eq. ( |2.18| ) , eq. ( p.23| ) does not define a geometric coordinate with respect 
to this derivative, meaning that d'^q^ ^ 5^ (actually, it is ill defined). However, g'^ is 
the geometric coordinate associated with the left derivative 9° with terminal Xq = 0, 



^ • °'"^ 9"gA* ■ r(l - a) Jo (xA* - x'^)" '^ 



(no sum over /i) 



By virtue of eqs. j^A^} and i ^Al\ ), 



W 



K 



(2.26b) 



(2.27) 



Since both derivatives ( p.26a| ) and ( |2.26bp will appear in the same theory (one in the 
action, the other in the equations of motion []TB|, section 2.3.5]), we are at liberty of 
choosing either in the action. In order to have a well-defined geometric coordinate 
system, the natural fractional differential is constructed via the left derivative ([|16], 
sections 2.4, 3.3] and references therein). 



(2.28) 
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with q^ and 9" given by eqs. ( p.23 ) and ( p.26b ). The symbol 9" may be actually 



regarded both as the partial fractional derivative with respect to x^ and as the one 



with respect to g'^, 9" = d/dg^. Therefore, we can define a geometric integral [jT6|, ^ 



r(a) Jo \x' J 



so that the action ( p.20 ) can be fully recast in geometric coordinate formalism: 

/+00 
d^qC. (2.30) 



2.1.3 Boundary 

Fractional spacetime Ai^ corresponds to the first orthant g^ > 0. Due to the 
presence of the boundary, translation and rotation invariance are globally broken, 
since these transformations would change the aspect of A^^ "looking from far away." 
However, as embedding transformations they do not affect the local symmetries of 
Ai^y which can still be investigated. 

This gives a sharp physical interpretation for the model defined by eqs. (|2.23| ), 



( p.26|) and ( |2.3CI| ). Despite being an integro-differential operator, the fractional 



derivative is an intrinsically local operator in the sense of fractional differential ge- 
ometry |6^, ^ (see also [0, section 3.1]). If A^^ possesses local symmetries and 
boundary effects are negligible well "inside" fractional spacetime, then the structures 
of f, (9" and d°' are mutually compatible. To show this, consider the one- dimensional 
interval [xo,a;i] and the behaviour of the differintegrals d^^,,. and 9"a.^, a 7^ 0, far 



from the terminals (see ||5^, sections 2.7.5 and 2.7.6, for left operators). The termi- 
nals are at x = xq and x = xi. For any finite xq, the asymptotic behaviour of d^^ ^ 
away from xq can be obtained by taking either x ^ 1 at fixed xq or —xq ^ 1 at 
fixed x: 

xo,x ~ \x, or ^xo,x ~ d^-xo,x- (2-31) 

Actually, the second asymptotic limit is tantamount to sending xq — )■ — cxd with x 
finite, and corresponds to the Liouville derivative. In our case, the lower terminal 
(or boundary) is finite and equal to x = xq = 0, and the first limit is the natural one 
for the fractional derivative inside the action. On the other hand, away from a finite 
Xi the right differintegrals behave as 

d:,x. '^' -"' 4% , or B:^^^ '^^1^"'^' d:^^^^^ . (2.32) 

For us, the upper terminal is x = xi = +00, and the second expression coincides with 
the Weyl differintegral in this limit. Thus, the theory defined by eqs. I\2.23i ), (^-20^ ) 



and (\2.3(\ ) can he physically interpreted as a local model of fractional spacetime where 
boundary effects are negligible. 
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What happens near the boundary? Consider now a D = 1 model where xi is 
finite: 



S'^=f dq{x)C. 
Jo 



Let f{x) be an analytic function on [xcXi] with integral singularities at x = xq and 
X = xi, i.e., such that it can be written equivalently as f{x) = {x — Xi)'^*fi{x), with 
/3j > —1, fi{xi) 7^ 0, and i = 0,1. Expanding f(x) in Taylor series around either 
terminal x = Xi, 



n! 

n=0 



and retaining the leading term, 






we obtain 



if a < /3o 

l)/o(xo) if a = /3o 
oo if a > /So 



lim (9"/)(x) = <^ r(« + l)/o(xo) if « = /3o , (2.33a) 



and 



if a < /3i 

lim (9"/)(x) = <( r(a + l)/i(xi) if a = /3i . (2.33b) 

oo if a > /3i 

Evaluation of differintegrals near terminal points lead to a collapse of their definition 
(even when the result is finite) and differential calculus thereon becomes inadequate. 
This conclusion is pleasantly expected in the light of the physical picture presented at 
the end of the paper. Roughly speaking, when all points of Ai^ are "too close to the 
boundary" the space and its boundary become one and the same. Such qualitative 
description reminds the topology of totally disconnected and post-critically finite 
fractals |B^, and in such regime fractional calculus, as a continuum approximation 



of these fractals, breaks down. Moreover, the use of Liouville/Weyl operators does 
not allow one to consider transient regimes because the limits |a;o,i| — ?■ cc correspond 
to looking far away from the terminal conditions, eqs. ( |2.31|) and ( p.32|) . To capture 



these regimes, one should take finite upper and lower boundaries. Later we will see 
that one can (and, actually, should) perform another extension of the model, in order 
to get phenomena which are transient in a different sense: promoting real fractional 
operators to complex ones. 
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2.1.4 Metric and distance 

Coordinate transformations between the Cartesian systems {x^} and a generic curvi- 
linear system {y^} are governed by the fractional generalization of vielbeins: 

ei := 5^g'(l/) , (2.34) 

which are D x D orthonormal matrices, e^e;^ = Sj. From these, one can define the 



fractional metric |65 



9,,u ■■= Vije^e^ , (2.35) 

where rju = diag(— 1, 1, ■ ■ ■ , 1) is the Minkowski metric. The fractional line element 
is then 

ds'' := [g^uidx^'y ® [dx^y]^ , (2.36) 



or, in geometric notation, 

ds^ = 9^^udq^ ® dq'' . (2.37) 

For ^A^, the metric is just the Minkowski metric, g^u = Vtiv The spatial distance 
between two points in fractional spacetime is then the 2a;-norm 



1 



A,(x,2/) := {[A(x^2/^)]"[A(a;^,y^)]"}^ := [Y^\x^ - y^\^- , (2.38) 



.M=l 



where A{x^,y'^) = |x^ — y'^\. This is a norm only if a > 1/2, i.e., when the triangle 
inequality holds. Therefore, we can further restrict a to lie in the range 



i<a<l 



(2.39) 



As already emphasized in |T^, one should not confuse eq. (|2.38|) with the choice of 



a p-norm (all topologically equivalent) in a given space: as a changes, so does the 
geometry of space. 

The last ingredient completing the definition of fractional spacetime is the Lapla- 
cian /C, entering the kinetic term of a scalar field theory living in A^^. Before dis- 
cussing this operator, we turn to two fundamental geometric properties of fractional 
spacetime: its fractal dimensions and the symmetry group. 

2.2 Hausdorff, spectral and walk dimensions 

In fractal geometry, there exist many definitions of "dimension" ||6^ (reviewed in 
[0). For ordinary spacetimes, there is no benefit in making distinctions among 
these definitions because they all agree in recognizing the topological dimension D 
as the number of degrees of freedom experienced by an observer in the measurement 
of distances, in diffusion processes, and so on. However, the topological dimension is a 
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bad indicator of the geometry of fractional spaces, and it is necessary to resort to the 
fractal machinery. This enters the picture in a simplified way, since the background 
is continuous. Full details are given in [|16]. 

The first useful indicator is the Hausdorff dimension ^h- Its operational definition 
on a smooth set is the scaling law for the volume V^^^ of a D-ball Bd of radius 5, 
which is 

g[BD{5)] = V^''\S)ocS''». (2.40) 

Then, in the limit of infinitely small radius, 

,„:.,JJ^. (2.41) 

5^0 In 5 

dn tells "how many directions" the observer feels in configuration space by making 



static measurements. From the scaling property (|2.25| ), one can already infer that 
dn = ao + {D — l)a < D. Fractional models are characterized by a Hausdorff dimen- 
sion smaller than or equal to the topological dimension D of embedding spacetime. 
A dynamical probe of dimensionality, on the other hand, consists in letting a 
test particle diffuse starting at point x and ending at point x'. In fiat fractional 
spacetime, this random walk is governed by a diffusion equation for the heat kernel 
P{x,x',a), 

(P^ - /Cf )P(x, x\ a)=0, P(x, x\ 0) = 6a{x, x') , (2.42) 

where a is diffusion time (a parameter not to be confused with physical or coordinate 
time), V^ is a diffusion differential operator of order f3, K,^ is the Laplacian (acting 
on the X dependence) defined on the Euclideanized ambient space, and 6a{x,x') = 
v~^{x)6{x — x') is the Dirac distribution in fractional space. The operators in ( p.42| ) 
can be chosen as 

V^^=d^,d^^,^d^,, (2.43) 

/C^ = 5^^ fd,d, + ^^^ , (2.44) 

with oq possibly different from the other charges. The kinetic operator is such that 
its eigenfunctions (Bessel functions of the first kind times a power) are the expansion 
basis of the fractional gemeralization of the "Fourier" transform |^. The heat 



kernel P{x, x, a) at coincident points x = x' and averaged in the space volume is 



called return probability. In fractional momentum space (see |]T6[ and section ^]5| 



V{a) oc a-^ . (2.45) 



The spectral dimension of fractional spacetime is then 



dmcr 

= (3[ao + {D - l)a] = (3dn . (2.47) 
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Finally, the walk dimension is the ratio (iw := '^d-a/ds- If di{ ^ ds, the diffusion law 
is said to be anomalous. For fractals, ds < d-a and d-w > 2, while diffusion with 
(iw < 2 is typically associated with jump processes. 



These matters were discussed in [16|, where the Hausdorff and spectral dimen- 
sions have been calculated for fractional isotropic Euclidean space. By definition, 
both the Hausdorff and spectral dimensions of a Lorentzian spacetime are calculated 
in Euclidean signature, so the results of [|T^ do not need any modification, except 
the separation of the fractional charge in time direction from the others. This poses, 
however, a caveat. Normal (or Gaussian) diffusion takes place if the derivative or- 
der of KP is twice the order of V^, i.e., when (3 = 1. So, Gaussian processes are 
produced by V^ = d^ and /C^ given above. Qualitative arguments suggest that 
Gaussian diffusion would be also achieved with V^ = d" and a kinetic term of the 
form such as /C^ = J2uoo^fi°'- ^^ f^i^ case, however, for ao j^ a diffusion would be 
anomalous, since anisotropy in the fractional charge induces non-trivial couplings in 
the Laplacian /C^ in order to match the scaling dimension of derivatives of different 
order. 

If A^^ is non- anomalous and /3 = 1, then the Hausdorff and spectral dimen- 
sions coincide, eq. ( p.47| ). In particular, d^ < du ii (3 < 1, while if /3 > 1 fractional 



Minkowski spacetime cannot be considered a fractal [|I6|. As a general result, on 
anisotropic models diffusion processes are anomalous. Note that the measure in mo- 
mentum space is the same distribution Qa as in configuration space, so the Hausdorff 
dimension of fractional momentum space is cin for any (3. 

2.3 Fractional Poincare transformations 

The scaling property of the Lebesgue-Stieltjes measure of fractional calculus, eq. 
( p.25|) , is associated both with a non-integer dimension and with a self-similar non- 



trivial structure. Self-similarity implies the same structure at all scales. So, fractional 
models should be compared with self-similar fractals, while models with dimensional 
flow are multi-fractal structures, which have different self-similar properties in ranges 
centered at different scales of magnification. However, in fractional spacetime the role 
of isometrics is far more significant than that of similarities, as discussed in section 
4.2 of |]IB|. Clearly, a particular presentation of fractional calculus breaks all Poincare 
symmetries, via the definition of measures which are neither translation nor rotation 
invariant. This may cause to believe that fractional systems, and fractal systems 
in general, are of little or no physical significance as classical and quantum field 
theories, where Lorentz invariance is an essential ingredient. For a generic model with 
absolutely continuous Lebesgue-Stieltjes measure dg{x) = v{x) d^x, it was shown 
that Poincare algebra was deformed even if the action itself is Poincare invariant 



IJ]. The key ingredient to obtain this result was to start with a Poincare-invariant 
Lebesgue-Stieltjes measure, i.e., regard v{x) as a Lorentz scalar. However, confusion 
arises when one chooses a particular profile for the measure weight v{x)\ any such 
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profile, which is needed as a concrete realization of anomalous scaling, explicitly 
breaks Poincare symmetries. For instance, a profile such as v{x) = |x|~'^'^~^''" is 
rotation invariant but neither translation nor boost invariant. Fractional models 



display explicit profiles v{x) by definition, and the attitude of [13|-115| is no longer 
tenable. Therefore, it is necessary to reexamine the issue of symmetries. 

In ordinary relativistic field theories, one observes that the line element (cross- 
product symbol ® omitted) 

ds^ = rj^^clx'^dx" (2.48) 

is preserved by the isometry group defined by the Poincare transformations 

x'^ = A^x" + a" , (2.49) 

where A(^ are D x D constant matrices such that 

a^,a; = A^,{A^,)-' = S^^ , 

and a'^ is a constant vector. Then, 

ds' = Tj^ydx'^dx''^ = 7] ^jjA^A^^dx^ dx" = rjp^dx^dx'^ . 

The Lorentz transformations A(^ = dx'^{x)/dx^ include spatial rotations and space- 
time boosts (det A = +1), as well as improper discrete transformations such as time 
reversal and parity (det A = —1). The A are actually frame transformations, since 
they act in internal space and vielbeins transform as e' = AjC^, but in homogeneous 
spacetimes one can choose one and the same frame for every point (the so-called Fermi 
frame), and tangent and Minkowski spaces are identified. 

One can construct actions which are Poincare invariant under proper transforma- 
tions. Lagrangians C[dx, ^'^{x)] of some fields </?* are defined to be proper scalars. On 
the other hand, the Lebesgue measure y/—gd^x is invariant, too, because det A = 1. 
In fact, in first-order formalism the volume element can be written as 

I det e\d^x = —6^° A . . . e^°~^ei^...i^_^ , (2.50) 

where e^ := e^ dx^ and 6^o--^d-i jg ^j^g Levi-Civita symbol, which is an internal pseudo- 
tensor: t'i^^...i^_^ = det AAj'^ . . . Aj°2l^.hy.JD~i- Applying a Lorentz transformation, 
one finds that \(\.eie'\d^x' = det A| det e|(i'^x. Thus, an action and equations of 
motions defined in a given coordinate frame {x} will be functionally the same in 
another frame {x'} related to the other by a proper Poincare transformation ( 2.49| ). 



Observers defined in local inertial frames experience natural phenomena governed by 
the same set of equations: physical laws depend neither on the position nor on the 
orientation of the observer's frame. 

In the case of fractional Minkowski space, the line element is given by eq. ( p.36| ) 
or, in geometric notation, 

ds^ = r^^^ dg" dg'^ . (2.51) 
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Consistently with eq. ( p.34 ), this suggests to define fractional Poincare transforma- 
tions which are hnear in geometric coordinates: 



q'^ = A^X + a' , 






(2.52) 



where a'^ is a constant vector and 



AM 

u 



d'^q'^ dq 



ii^ 



Qaqu Qqu 



(2.53) 



are D x D constant matrices, associated with the ordinary SO{D — 1,1) group. 
Equation ( p.52[ ) is in accordance with [^ section 4.2], where it was argued that 
fractional Euclidean space Sj^ is characterized by the group of affine transformations 



g'^ = A^g^ + fi^ . 



The line element ( 2.51|) is preserved under the transformations ( |2.52|) . Crucial 
to this result is the fact that the fractional differential, made of Caputo derivatives, 
is zero on a constant. If we had used Riemann-Liouville calculus, we would not have 
been able to write ( 2.52| ) in such a simple form. 

We now discuss whether eq. ( p. 521 ) is the only transformation preserving the 
fractional line element for general q. Let us first recall a textbook exercise for the 
integer case |^ . Under a general non-singular coordinate transformation x^ — )■ x' 



'M 



the line element (|2.48| ) changes as 



ds' = rj^^dx'^dx'" = ri^^dpx'^dax"'dx''dx'^ 



„/2 



but imposing ds' = ds , one gets i^pa = rj^ydpX dax' . Differentiating with respect 



to x^, one has 







V^lu 



2^11^ 



d^x'^ dx'" dx'" d^x 
+ 



dx'^dxP dx'^ dxP dx'^dx" 



adding and subtracting the same equation with, respectively, t -H- p and t ^ a, one 
obtains an equation with six terms. They cancel one another except two identical, 
giving twice the first term. Since rj^^, and dx'^ /dx" are non-singular, one ends with 

the condition 

Q2^n^ 



0, 



dx'^dxP 

whose solution is eq. ( 12.49] ). In the fractional case, the same argument does not go 
through as smoothly. The failure of the Leibniz rule 



df 
dx[f{g)] = ^^^9 



(2.54) 
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and its replacement with 

^ ^^^ ^ r(l-a) (x-xo)° 

+ g sm[^(,-.)]r(l + .) _ ,^^.-.(,.^)(^) (2.55) 

^ 7r(j-a) r(l+j) 

forbid mixed derivatives to combine in a simple way. At most, one can recognize 
eq. (|2.52|) as a sufficient but not necessary condition for the line element to be in- 
variant. The fact that we have infinite terms all of different order, however, makes 
it likely that eq. ( |2.52| ) is also necessary, unless miraculous cancellations take place. 



Just like for integer transformations, the orthogonality relation in ( p^.52D implies 
that 



D-l 



CKY = 1 + E(^°)' ^1' (detA)^ = 1, (2.56) 

SO one can distinguish between proper and improper transformations. We dub the 
semi-direct product of translations and fractional Lorentz transformations on the 
embedding coordinates x'^ the fractional Poincare group IIq,. 

Clearly, fractional transformations A(^ neither act linearly on embedding coordi- 
nates x^ nor are simply given by the elements of A^ to the power of a: (rfx"^)" = 
(A(^(ix'^)" under an integer transformation, while dq"^ = A^dg*^ = A(^((ia;'^)". We can 
find an approximate relation between integer and fractional Lorentz transformations 
when a = l — e, e^l. Noting that 

q = x + x{l-j-\nx)e + 0{e'^), (2.57) 

where 7 is Euler's constant, fractional Poincare transformations reduce to integer 
ones up to correction terms: 

a;; = [1 + (7 - l)e] a;; + 0(e^ x in x) . (2.58) 

The fact that fractional symmetries are intrinsically non-linear resembles the situa- 
tion in non-commutative /t-Minkowski spacetime, which is invariant under deformed, 
non- linear Poincare symmetries |]68[-|7^. The relation between fractional and non- 



commutative theories is discussed elsewhere 

Because of the boundary at x^ = 0, the fractional Poincare transformations 
( p.52|) are not global symmetries of fractional Minkowski spacetime. Yet, they are 



the guiding principle to write an invariant action up to boundary terms. We show 
first that the fractional measure, derivative and differential are invariant under the 
fractional isometry group Ua, up to boundary terms. The task is easy in geometric 
coordinates. The proof for the integration in eq. ( |2.20| ) is formally the same as 
in ordinary calculus, with e^. A/ and x in eq. ( p^.50| ) and below replaced by their 
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fractional counterparts. The lower terminal in ( p.20 ) is translated by a^ but this 



does not affect the physics well inside the spacetime bulk, according to section |2Tl . 
For consistency, writing the measure as in eq. ( |2.30| ) should yield the same result. 
To check this, it is sufficient to verify that the fractional Caputo derivative transforms 
as a covariant vector. This is true, thanks to eq. ( |2.53| ): 



' Kl^,- (2-59) 



In particular, operators of the form rj^^ (pd'^d'^ cf) are Lorentz covariant in a fractional 
sense, if is a scalar. 

Looking at the definition of 5" this result would not have been obvious. For 
< a < 1, the left derivative can be written in terms of the geometric coordinate q 
(index fi omitted): 

1 r dr' 

(dV){x) = — -/ ^^— 9.,/(x') 

r(l - a) Jq {x - x'Y 



'^ r / — dyf(x — y) 

r(i - a) J, r ' ^ ' 

YTT^ r%9,,f{q,q'), (2.60) 



where f{q,q') = /{[r(l + a)qY^"' — [r(l + a)q'Y^"'}. For the Weyl derivative, one 
obtains a similar expression with the upper terminal x replaced by +oo and f{q, q') = 
/{[r(l + a)qY'°' + [r(l + a)q'Y'°'}. Performing a fractional Poincare transformation 
and changing integration variable leads to a seemingly no n- vectorial expression, just 
like acting naively on d/dx^ with eq. ( p. 491 ). 

To summarize, the explicit coordinate dependence x^ of fractional operators is 
just a presentation in the embedding and fractional covariance must be defined in 
the space of geometric coordinates q^. 

3. Dimensional flow 

So far, we have discussed various ingredients for the construction of a spacetime 
characterized by a measure of fixed fractional order. This was done for the purpose 
of simplifying the description of an unfamiliar type of geometry to the bone. However, 
a major goal of this proposal is to introduce the often-advertized idea that geometry 
changes with the scale or, in more colorful words, to give spacetime a multi-fractal 



structure. In ||T3[-||T5[, it was assumed that the parameter a would somehow flow 
from some finite a = a*<ltoQ; = l, without however giving any detail about how 
this flow takes place. 

One could continue to keep a fixed and develop the concrete example of a scalar 
theory, were it not for a simple but maybe surprising fact. Namely, a multi-fractional 
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classical structure is intimately related with the quantum structure of this class of field 
theories, and it determines its renormalization group properties. The reason is that 
the presence of all possible fractional operators in the classical action determines a 
hierarchy of scales. This hierarchy can be interpreted both as the very definition 
of multi-fractality in the fractional context and as a self-consistency requirement for 
renormalization. Not only these interpretations are not mutually exclusive, but they 
practically amount to one and the same. Hence the logical necessity to first present 
the definition of mult i-f r act ionality, and then construct consistent classical actions. 
Here we begin to carry out this programme. 

One possibility would be to promote a to a Lorentz scalar, a function of spacetime 
coordinates: 

a -^ a{x^) . (3.1) 

The dimension of spacetime would become dynamical, a possibility already consid- 
ered in the past ^^. One should add a kinetic term and a potential for this new field 



into the total action, which would be augmented by an extra contribution 5'[a(a;)]. 
However, the system would quickly become intractable due to the complicated in- 
tegration measure, not to mention the transcendental dependence on a in ^[^(a;)]. 
We also ask ourselves whether we wish a model to become fractional at small spatial 
scales or at early times, or both: the first case is a universe becoming fractal be- 
low a certain critical scale at any given time, while the second corresponds, roughly 
speaking, to a spacetime becoming fractal near the big bang. Here we face physi- 
cally inequivalent scenarios, depending on how we define the action. For example, a 
realization of a spacetime becoming fractional at early times would be to allow the 
fractional order of time direction to be ao = 1, while that of spatial directions to 
be only time dependent, a^ — )■ a{t). Then, one would have no problems in integrat- 
ing, if not for the fact that time and fractional spatial derivatives would no longer 
commute. If only ordinary derivatives were present in the classical action, the only 
minor complication would be to integrate the measure by parts when deriving the 
equations of motion. 

3.1 External scale or renormalization group picture 

The main problem arising with eq. ( p.l| ) is that geometric coordinates would be 
defined through transcendental expressions and, in general, the fate of fractional 
Poincare transformations is not clear. In a Lorentz-covariant framework, the presence 
of a characteristic scale can spoil the symmetries of the system, unless one introduces 
the latter with care. This suggests that treating a as a field may be unsatisfactory. In 
alternative, one can adopt the perspective of critical systems and regard a as an order 
parameter. More precisely, one can parametrize a not with spacetime coordinates, 
but with an external scale parameter governing dimensional flow. Physically, one 
can introduce a critical length/time scale i^ below which the system flows to a deep 
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fractional regime, and above which it occupies the whole embedding space. Assuming 
that a acquires a finite critical value a^ at the bottom of this regime, we have 

"^^) ^ 1 1 as £ » 4 • ('•') 

Notice that a similar behaviour would appear in a renormalization group picture, 
where the parameter a runs with the energy scale. If interpreted as a fundamental 
scale (which is true in fractional theories of real order), £* may be associated with 
the Planck scale. In this context, as in modern approaches such as asymptotic safety 
gravity, it does make sense to consider scales below £*, because the latter is not a 
cut-off of the theory. Anyway, we refrain from the identification l^, = ipi because in 
the complex-order theory the scale hierarchy will go through a little revolution, and 
the Planck scale will be pushed further deep in the UV spacetime structure. 

Defining the adimensional multi- fractional external time a := i/i:^, we can re- 
place the limits in eq. ( |3.2| ) with 

. N fa* as or ^ 

aia) ~ < , . (3.3) 

^ ^ [1 as a ^ +00 ^ ^ 

For instance, profiles realizing ( p.3| ) are 

a(cr) = — ^ , q;(o") = a* +(1— a*) tanh cr , a(cr) = «*+(!— a*)erfcr , (3.4) 

1 + CT 

where erf is the error function. The latter is solution of the diffusion equation in 
Euclidean space with the Heaviside distribution as initial condition. 

Let eq. (p.20|) be an action characterized by a measure ga{x) of dimension di{{a). 
In particular, dii{a) is given by eq. ( p.47|) . In the external time picture, the funda- 



mental action is decorated with an extra integration over external time, 

(3.5) 




where g{a) is a one-parameter measure. Since a is an external time parameter, 
spacetime covariance is respected if Sa is covariant. The spacetime whose measure 
is a superposition of fractional Minkowski measures will be called multi-fractional 
Minkowski spacetime, denoted by A^f . 

3.2 Multi- fractional spacetime 

The external-scale Ansatz ( p.5|) admits a neat interpretation in terms of multi-fractal 
geometry. Just like fractals, multi- fractals do not have a precise definition but, in- 
tuitively, they are sets with scale-dependent fractal properties, on which mass dis- 
tributions (i.e., measures) do not obey a simple power law (see, e.g., [BO, chapter 
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17] and [^). We have seen that the Hausdorff dimension of a smooth set is con- 
veniently determined by the scahng law of the measure of balls of infinitely small 
radius, eq. ( |2.41| ).^ One can make this definition local and dependent on the center 
X of the ball, provided the following limit exists: 

di^ := hm -— , (3.6) 

5^0 In 5 

a quantity called local dimension or Holder exponent of g at x. Consider now the set 
J-rf of the points where the local dimension exists and equals d: 

J^^=[xeMP : di^ = d] . (3.7) 

In other words, ball measures centered at every x in Fd scale as V^^^((5) = q[Bd{x, 5)] ~ 
5*^ as 5 — >■ 0. As (i varies, so does the set Fd and its Hausdorff dimension. The prob- 
lem of calculating this dimension, called fine multi-fractal (or singularity) spectrum, 
is the subject of fine multi-fractal analysis: 

h{d) := rfH(^d) . (3.8) 

The fine spectrum /h gives information about the limiting behaviour of q at any 
point. 

Self-similar sets well illustrate the methods of multi-fractal spectral analysis. We 



recall some basic definitions from [^ , discussed in |]T6[ in the context of fractional 
spaces. Consider a set of A^ maps Si : MP — > R^, 'i = l,...,A^>2, such that 



^[Si{x),Si{y)]<\A{x,y), x,i/eM^, < A^ < 1 , (3.9) 

where the distance A between two points is A{x,y) = \x — y\ in ordinary integer 
geometry. Any such map is called contraction and the number Aj is its ratio. If 
equality holds. Si is a contracting similarity or simply a similarity. By definition, 
self-similar deterministic fractals are invariant under contraction maps and can be 
expressed as the union of the images of Si []7B| , 



N 



T=[jS,{T). (3.10) 



i=l 



Suppose the strong separation condition holds, i.e., there exists a closed set U such 
that Si{U) C U ior all i = 1, . . . , N and Si{U) U Sj^i{U) =0. J" C f/ is constructed 
taking sequences of similarities and the intersection of sets Uk = Si-^ o ■ ■ ■ o Sii^{U). If 
\U\ = 1, then the diameter of the k-th iteration set is the product of similarity ratios, 
l^fcl = Ajj . . . Aj^. Let < Qi < 1 he N probabilities (or mass ratios, or weights) such 



^For a non-smooth set embedded in M-^, eq. ( p.4lD still defines a "fractal dimension" but, in 
general, it does not coincide with dfj- 
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that Yli9i ~ -'-• O'^^ "^^^ imagine to distribute a mass on sets Uk by dividing it 
repeatedly in A^ subsets of Uk, in the ratios gi : ■ ■ ■ : g^. This defines a self-similar 



measure g with support J^, such that Q{Uk) = gi^ . . . gt^. and, for all sets A C J^ [[f^ , 

f,(A) = 5^^,f,[5ri(A)]. (3.11) 

The case A^ = +00 defines so-called infinite self-similar measures, describing fractals 
with an infinite number of similarities ^^ ^. Given a real parameter u, we define 



the singularity (or correlation) exponent 9{u) as the real number such that [[79|-p2| 

N 



^5rA?") = l. (3.12) 



The correlation exponent exists and is unique, since < Aj, f^j < 1. As a function of 
u, 6 is decreasing and limu_j.±oo^(^) = =Foo. Finally, the generalized dimensions are 

d(u):=^^, M^l, (3.13) 

M — 1 

and a non-singular definition, which we do not report here, is employed for u = 1. 
Fractals characterized by just one dimension at all scales are special cases of multi- 
fractals. For deterministic fractals, the probabilities are all equal to gi = 1/N. In all 
fractals with equal contracting ratios Aj = A, the generalized dimensions all coincide 
with the capacity dc of the set (e.g., [0), which is also the Hausdorff dimension. 
From eq. (|3:i2|) , 

A^W In A^ 

iV^^ = l ^ diu) = —^=:dc = dn. (3.14) 



For self-similar measures ( p.ll| ), and under some weak assumptions, the spectrum 



enjoys a number of properties |BB]. First, the support of the fine spectrum is in 
a certain finite interval [(imin, c^max], on which /h is given by the Legendre transform 
of ^: 

M{d) = mi[e{u) + du]. (3.15) 

For d ^ [rfmin, c?max], J^d is the empty set. Second, /h is a concave function of d. The 
maximum of the spectrum is at d = d{0) and equals the dimension of the support 
of the measure, maxaf-aid) = f-ii[d{0)] = 6{0) = dnisuppg). In this case, eq. ( p.l2| ) 
becomes the well-known algebraic condition for the Hausdorff dimension of a self- 
similar set: 

N 

^ ^dH(supp,) _ ^ _ (3.16) 

The support of the measure can be a fractal; for instance, the multi-fractal Cantor set 
is a mass distribution over the Cantor set, which is its support. Third, at m = 1, for 
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eq. ( p.l2|) ^(1) = and the spectrum equals the Hausdorff dimension of the measure: 

fn[d{l)]=d{l) = dn{g).^ 

There are a plethora of physical phenomena described by multi-fractal geometry. 
Some multi-fractal systems can be approximated by fractional dynamics. For exam- 
ple, in the case of the fractional Fokker-Planck-Kolmogorov equation, one simply 
replaces fractional derivatives with sums of derivatives of any allowed fractional order 
]. Multi- fractional Brownian motion is another instance [^, |7^, |5^-||87||. 

Let us go back to fractional spacetime Ai^ of fixed order a, and consider for sim- 



plicity the isotropic case a^ = a. In [|I6|, we have argued that its Euclidean analogue 
S^ can be characterized by two similarity maps acting on geometric coordinates. For 
any direction /x, 

Si(g^) := Xq" , S2(g^) := (1 - \)q^ + A , (3.17) 

where A is arbitrary and chosen to be the same along all directions. Since A is 
arbitrary, this set is trivially self-similar; yet, indeed it is self-similar. A mass would 
be equally distributed on, say, N = 2^ subsets, with probabilities Qi = 1/N. In 
geometric coordinates, the scaling is A = 1/2 = N'^^/^ . Applying eq. ( p.l4|) , c^h = 
— In N/ In A = -D, and c^h = Da for embedding coordinates. From the point of view 
of the fractional embedding spanned by the x, the scaling ratio is 

and the probabilities can be regarded as a-dependent for a given A: 

(?. =^«:=A^". (3.19) 

The same conclusion holds for fractional Minkowski spacetime. 

Collecting these results, we determine the extension of the measure Qa to a 
multi-fractional measure with the self-similar structure ( p.ll ). A non-degenerate set 



of probabihties {gi} must be introduced to weight the sub-copies of A^^ differently, 
with the index i running on a given set (with more than two elements, possibly). 
Each copy is labelled by a, which plays the role of generalized dimension at a given 
scale. Thus, and by eq. ( p.l9| ), the labelling of the probabilities is assumed by a. For 



discrete a, one should exchange the a-order integration in the fractional model with 
a sum of integrals over all possible ranges in a: 

(3.20) 



S-- 




^a } 


1 
2 


< 


a 


< 1. 



*^The Hausdorff dimension of a measure is 

dniQ) '■— ini{diiE : E is a, Borel set with g{E) > 0} . 

Intuitively, it is the dimension of a set on which a significant part of the mass distribution g is 
concentrated Ipffl. 
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The external time/RG picture ( |3.5|) is nothing but a continuum version of ( p.20| ). 

The coefficients Qa are probabihties from the point of view of fractal geometry, 
and coupling constants from the perspective of field theory. Assuming the Lagrangian 
to be the same for all contributions, the g^ are dimensionful in order for S to be 
dimensionless. Then, they determine the scale at which geometry changes. Consider 
a simplified isotropic model with one such (length) scale £*, where D-dimensional 
integrations (whatever the choice of fractional calculus) are given by 

/D = /S^+^f"^-"^^/S% [lD] = -Dau i<ai<a2<l. (3.21) 

The measure is binomial and volumes are made of two pieces. The volume of a D-ball 
of radius i? in a space endowed with this structure is 



V^\R) = i 



Dai 



T3 \ Da\ / f? ^ ^'^'^ 



^D.a. ^ + fi 



D,ai I "^ "T i'D,a2 



(3.22) 



where flD,a is the volume of a unit ball. Depending on the size of the ball (i.e., on 
the scale one is probing), either term will dominate over the other, thus giving two 
types of scaling laws. For a small ball {R < Q, V^^^ ~ R^"\ while V^^^ ~ ^^"^ 
for a large ball (i? :» £*), where R = Ri^ 01/02 ^^ ^j^^ radius of the ball measured in 
macroscopic length units (units effectively change with the scale in a multi-fractional 
setting). 

3.3 Dimensionality and the role of dn = 2 and D = A 

An important consequence of the geometric picture outlined in section p.l.4| is that, 
for a given topological dimension D > 1, not all fractional measures are possible. In 
the (at least) spatially isotropic case, the Hausdorff dimension d}i = Oq + {D — l)a 
is associated with 2a-norms only if a, ao > 1/2. If ao 7^ 1, this implies 

D <2dn, (3.23) 

while for ag = 1 one has 

D <2du-l. (3.24) 

For instance, dimensions di{{a) < 2 are not associated with 2Q;-norms, unless D < 3. 
The interesting dimension d}i{a) = 2 is achieved at the critical value 

Imposing a* > 1/2, the critical value a* exists for D < 4 in the fully isotropic case 
or in the spatially isotropic case with a^ ^ 1. This excludes an integer time direction 
if a two-dimensional regime is required. See table |I|. 
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Furthermore, if one also assumes that a, ao < 1, there is a lower bound for D, 



D>dn 



(3.26) 



stating that fractional spacetime is embedded in the abstract space M^. Thus, 
the dimensionality of the critical point with the requirement of the existence of a 
natural norm provides a guiding principle in determining the maximal dimension 
in the infrared. A less clear indication of this feature was noticed in |14] for the 
structure of a scalar propagator. The reader may wonder whether one can also 
produce geometries where D < du, a^s it happens in certain quantum gravity models. 



In section 3.3.3 we will comment on this. 





dn<2 


dn = 2 


du>2 


ttO < 1 
«o = 1 


D<3 
D<2 


D <A 
D<3 


D < 2dn 
D<2dn-1 



Table 1: Allowed topological dimension of the embedding for a fractional spacetime with 
Hausdorff dimension dii{ao,a) and natural 2Qi-norm. 



In the embedding picture considered in this paper, one can now reinterpret the 
parameter a also as a measure of maximal dimensional dispersion. In fact, the 
difference between the maximum and minimum topological dimension for a given du 
is ADmax = 2(iH — c^H = c^H = Da for the fully isotropic model, so that 



AD, 



D 



a . 



(3.27) 



On the other hand, the inverse of a is the dispersion between the maximum and 
minimum allowed Hausdorff dimension, Adu = D — = D = d}i/a: 



Adn 
dH 



1 
a 



(3.28) 



Also, the difference between the maximum allowed topological dimension and the 
actual D is AD = 2du- D = D{2a - 1), hence 



AD 
~1T 



2a- I. 



(3.29) 



The special role of D = 4 consists in allowing, in the same dimensional flow, a phe- 
nomenologically viable macroscopic scenario and an ultraviolet configuration with 
natural 2a-norm and Hausdorff dimension equal to 2. The requirement of having 
a geometric norm throughout the dimensional evolution is crucial. One should be 
careful in drawing the conclusion that D = 4 is favoured, since d-a = 2 was only 
suggested by quantum gravity arguments (and by a preliminary analysis of the UV 
finiteness of the theory in D = 4, section |4.(j|). However, there exists a mysterious 
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relation between the fundamental constants of Nature which further supports the 
very special role oi dn = 2 as the characteristic dimension of a quantum theory of 
gravity. Including Planck's constant h, the electron charge e, Newton's constant G, 
and the speed of light c, one can construct a dimensionless constant in D dimensions 
as C = h?~^e^^^G^^^c^^^ |Q. The dimensional arguments leading to this combi- 



nation are unchanged if one replaces the topological dimension with the Hausdorff 
dimension. Replacing also h with the Planck length 



'PI 



no 



c 



3 



1.6163 X lO^^^m, (3.30) 



the same constant can be recast as 

C = 4f-"«)e'^H-2G^-ic2(2-'^«) . (3.31) 

In (in > 3, the Planck length appears in negative powers, a hint that perturbative 
quantum gravity is non-renormalizable in these dimensions. At (in = 3 the Planck 
length disappears (and, in fact, perturbative gravity is renormalizable in three di- 
mensions). Remarkably, in (in = 2 the fundamental constant coincides with (the 
square of) the Planck length, C = ip^, while all the other couplings disappear. 

This argument highlights the peculiar status of c^h = 2 in quantum theories of 
matter and gravity.^ In the present context of dimensional flow, we can take it as 
a hint of the special role of D = 4. Other observations select D = A topological 
dimensions as special |^. Just to mention a few examples (which all assume one 
time direction), stable planetary orbits appear only in D = 4; stable Bohr atoms, 
in D < 4; in D = 4 the number of generators of spatial translations and rotations 
is the same, with consequences for electromagnetism; wave signals propagate free of 
reverberation and distortion only in D = 4; spacetimes with D < 3 do not contain 
gravitational waves; chaos may appear in dynamical systems only in D > 4; and 
others. It would be interesting to embed these arguments in a multi-fractional model 
and apply them to the Hausdorff dimension. 

3.4 Experimental bounds on dimensional flow 

The Hausdorff dimension is a physical observable. The volume law of, say, a mass 
distribution concentrated in a spherical region of size ~ i? is represented in a log-log 
plot, 

\nV{R)=\nnD,a + Da\nR, (3.32) 

where flD,a is the volume of a D-hall with unit radius. One can, in principle, obtain 
independent measurements of volumes and radii and determine both the dimension- 
ality of the mass distribution (tilt of the line) and the solid-angle factor (offset). With 

''Noting, as it is usually done, that [G] ^ 2 — D vanishes in two dimensions and then deducing 



that RG properties of gravity are special in I? = 2 is not quite the same. Beside gravity, eq. (3.31) 
involves also other forces of Nature. 
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adequate technology, one can translate this type of experiments in the realm of space- 
time geometry. Experimental constraints on fractional models can be obtained from 
observations where gravitational effects are almost or completely negligible. These 
can be, mainly, of three types: particle physics tests, equivalence principle/Lorentz 
invariance tests, and post-Newtonian tests. We saw in [|1^, section 3.5] that different 
presentations of fractional geometry yield different output values for, e.g., volume 
measurements, but they enjoy the same scaling properties for a given a. Thus, if we 
assume that at scales about and above those probed by accelerators we are already 
out of the multi- fractional regime, we can consider a d}^ = D — e expansion for a 
fixed-a configuration, and focus our attention to the parameter e. This falls short of 
constraining the characteristic scale £* at which the UV critical point is approached, 
since these experiments correspond to scales £ ^ £*. Anyway, later we will also give 
a bound on i^, from particle physics. 

In table ^we report the expressions of the unit volume flD,a at D = 2, 3, 4, where 



a 



'-D- 



(3.33) 



Calculations where du = D — e correspond, in fractal geometry, to regimes with low 
lacunarity, i.e., where fractal space is almost translation invariant 



D 


^D,l~e/D 


^D-e,l 


2 
3 

4 


7r(l-0.42e) 
f (l-0.54e) 
^(l-0.63e) 


7r(l-0.36e) 
f (l-0.22e) 



Table 2: Volume 0,£) „ of unit D-balls in various dimensions, for a ~ 1. The corrections 



in traditional dimensional regularization are shown in the last column [16|. 



The coefficients in front of e are all of the same order of magnitude in fractional 
models and in spacetimes with dimension modified according to the dimensional- 



regularization scheme ||90|-||92[, so we can accept bounds on the latter and apply 
them to the former. For instance, measurements of the anomalous magnetic mo- 
ment g — 2 oi the muon can compare the theoretical prediction in four dimensions 
with that in spacetimes with dimension modified according to the dimensional- 
regularization scheme. The order of magnitude of the upper bound on e was es- 
timated as |e| < 10~^ in [^, while from the anomalous magnetic moment of the 
electron, e ~ lO^l^fthcc 



fi'cxpl Uni m and one finds |e| < 10" 



experimental determination of the muon g — 2 |9^, where gcxp — fi'th. 
have 

lei < 10"^ £-10"^^ 



9^ . Taking the latest 
^ 10-11 



we 



m. 



(3.34) 

Going to atomic scales, experimental measurements of the Lamb shift in hydrogen 
yield p], || 

(3.35) 



< 10 



-11 



10^11 m. 
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tighter than the previous bounds. 

As far as intermediate-scale experiments are concerned {i < 10^ m), one could 



reinterpret tests of the equivalence principle ^7j. Of course, at mesoscopic scales 
above particle physics and below planetary we know that the equivalence principle 
holds with great accuracy, and local physics is described by Minkowski spacetime 
(plus eventual corrections of general relativity, measurable already at the size of the 
LHC accelerator). However, it is a legitimate question to ask what the bounds would 
be on a 4 — e geometry. So, one can use fractional equations correcting Euclidean 
geometry and apply them to the data of these experiments. The equivalence principle 
is closely related to Lorentz invariance, and tests of the latter would provide parallel 



constraints on fractional effects p8| , |99| . Upper bounds on e, and hence on the 
dimensional flow at such mesoscopic scales, should be quite stringent. 

These bounds should be compared with others done at larger scales. Anomalous 
correlation functions result in deviations from Newton's law and a precession of 
planetary orbits. Taking into account the theoretical prediction of general relativity 
for the precession of Mercury, any dimensional effect should be smaller than the 
experimental error, thus yielding |^ ^ [LOCI|] 



|e| < 10"% f-lO^^m, (3.36) 

less tight than eq. (|3.35|) . A similar bound is obtained from pulsar measurements at 
a distance 6 ~ 10"^ ly from us |p,OCI|| , so that one can regard ( |3.35| ) as valid not only 
here and today, but also in a spacetime sphere of size 6. 

Finally, a fit of the black-body spectrum of the cosmic microwave background 
gives an upper bound on e at spacetime scales comprised between the decoupling era 
and today ||101|| : 

|e| < 10"^ £~14.4Gpc. (3.37) 

The best-fit value of e is strictly positive, in agreement with the direction of dimen- 
sional flow in fractional theories. 

All the above estimates rely on a number of assumptions which require a careful 
scrutiny in the present framework of fractional spacetime. Among these assump- 
tions, we mention naive implementations of dimension effects (via dimensional reg- 
ularization, which is of a non-dynamical nature), the use of unmodified Einstein or 
Schrodinger equations, and integer time direction (in our language, fractional charge 
ao = 1). If the parameter e is rendered dynamical, as it would naturally be in multi- 
fractional dimensional flow, then the above particle-physics bounds are no longer 
reliable and, in fact, a determination of spacetime dimension becomes much more 
difficult. Examples are measurements of oscillations of neutral B mesons and of the 
muon g — 2 ||102|| : At mass scales M > 300 -=-400 GeV, any 2 < du < 5 is compatible 



with experiments. This translates into a rough upper bound for £^: 

4<10"i®m. (3.38) 
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To carry out a complete revision of the results presented in this section, one will 
have first to extend the model to curved spacetimes with gravity. This is part of the 
theoretical programme here proposed for future studies; in parallel, the improvement 
of data could allow an update on some of the above results. Lamb shift measurements 
seem to be the most promising for the tightest absolute constraint. However, one 
should not underestimate the bound ( |3.37D , which, if confirmed, would constrain the 
end of the multi-fractional era of the Universe at times prior to the formation of the 
cosmic microwave background. 

4. Scalar field theory 

A Lorentz scalar is a favourite guinea pig with which to check the properties of a 
field theory in a given geometry setting. We studied a real scalar field in WM for 



a generic absolutely continuous Lebesgue-Stieltjes measure. Here we do just that 
but focussing on fractional measures and power-counting renormalization properties. 
The latter are important when fractional theories are regarded as fundamental, in 
which case one should check that they are ultraviolet finite. 

4.1 Power-counting renormalizability 

A standard power-counting argument is sufficient to understand qualitatively the re- 
lation between coupling dimensionality and renormalization properties of a field the- 
ory ||103|| - |T05| . The reader already acquainted with it can skip this section. When 



constructing perturbation theory, one must take into account all possible gauge- 
inequivalent interactions order by order in the effective low-energy action. Some of 
the couplings diverge when the regulator in the regularization scheme is removed. 
However, if the operators O ~ J dgOd associated with these couplings (/ren are al- 
ready present at the tree level, with bare couplings g of dimension d^ — d (where 
c? > is the scaling dimension of Od)-, one can absorb the divergence into an effective 
coupling which is defined to be finite when the regulator is removed. Contrary to 
the bare couplings g, the effective couplings g^^ = g + gren are what one physically 
measures. If this procedure works order by order, the theory is said to be perturba- 
tively renormalizable, and hence physically predictive. This means that the number 
of physical couplings we measure at any perturbative order is finite. 

In the renormalization group picture, the physical action stems from the bare 
action by integrating out momentum modes greater than a certain energy cut-off 
scale E, and then removing the cut-off. An operator is said to be relevant if its 
associated coupling g has positive scaling dimension. On the other hand, operators 
with dimensionless coupling are called marginal, while operators with [g] < are 
irrelevant. In terms of dimensionless constants 

~g = gE''-''-, [g] = 0, (4.1) 
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the operator O associated with g scales as 



^ \ d-dH 



O = gJd0O,r^gi-\ , (4.2) 

where k is the momentum. Therefore, relevant operators are important at low en- 
ergies {k/E -C 1). Marginal operators are equally important at all scales. Their 
detailed behaviour is not obvious and these can be, case by case, either marginally 
relevant or irrelevant. Irrelevant operators become important in the ultraviolet 
{k/E ^ 1) but, contrary to what the name suggests, some of them can also alter 
macroscopic physics. Since, typically, there is a finite number of relevant opera- 
tors and also of marginal operators, macroscopic physics is described only by few 
observables. In fractional theories, this would not be the case if we included the 
infinite class of operators with fractional derivatives. Imposition of symmetries will 
drastically reduce this infinite multiplicity. 

If divergences are present, they correspond to local operators of dimensionality 
increasing with the order of the perturbation expansion. Suppose the bare action 
S contains only relevant operators; then, only a finite number of relevant operators 
(those which did not appear in S) will enter the effective action, and any divergence 
will be absorbed in the finite number of couplings {g}. For instance, in electro- 
magnetism the electron mass and charge have non-negative dimension in natural 
units, and one can formally absorb the divergences just in these two coupling con- 
stants, which are then determined by experiments. Conversely, if even one irrelevant 
operator appears in S, one can construct new irrelevant operators at each order. 
Explicit calculations can determine whether their couplings are finite or not. If they 
diverge, the perturbative approach looses predictivity because we can absorb all the 
divergences only by adding an infinite number of operators to the action. 

A theory is said to be power- counting renormalizable if 

[g] > (4.3) 

for all bare couplings g. This condition is not sufficient to guarantee that the theory 
be renormalizable in the sense of the full renormalization group flow, but it provides 
a good guiding principle in many situations. If a model is not power-counting renor- 
malizable, then it will likely be non-renormalizable unless remarkable divergence 
cancellations happen. An example is ordinary and supersymmetric gravity, where 
these cancellations do happen ||106|| and explicit calculations are necessary to settle 



the issue. 

The relation between the good UV behaviour of a theory and the absence of irrel- 
evant operators can be understood by looking at the superficial degree of divergence 
of a Feynman diagram. Consider a one-particle-irreducible FejTiman sub-graph with 
L loops, I internal propagators and V vertices. The superficial degree of divergence 
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6 is the canonical dimension of all these contributions: given a UV energy cut-off E, 
the divergent part of the diagram scales as E^. li 6 = 0, one has at most superfi- 
cial logarithmic divergences and the theory is power-counting renormalizable. When 
S < for every sub-diagram in a Feynman graph, the graph is convergent; if only a fi- 
nite number of Feynman diagrams diverge superficially, the theory is power-counting 
super- renormalizable. 

We can count divergences in the case of an ordinary scalar field theory in D 
dimensions: 

S= I d^'x [-09^9^0 - -m'(f - ^(f\ . (4.4) 

Each loop integral over momenta gives [d^k\ = D, while the propagator G{k'^) = 
— l/{k'^ + rn^) has [G] = —2. Interaction vertices do not carry dimensionality and, 
overall, 6 = DL — 21. Since I > L, the maximum superficial degree of divergence 
can be L{D — 2). L is the number of independent momenta, given by / minus the 
number of relations they satisfy among themselves: these are V — 1 (one for each 
vertex, up to the total momentum conservation), so that 

L = I-V+l. (4.5) 

This result is often called Euler's theorem for graphs. With only mass and a 0^ 

interaction, for each vertex there are p lines, so that pV = N + 21, where A^ is the 

number of external legs in the diagram. Replacing L and / with these expressions, 

one obtains 

D 

~2 



S = DL-2I = D- [g,]V - ( - - 1 ) iV, (4.6) 



where we used [gp] = D—p{D—2)/2. This formula can be also derived by dimensional 
arguments. A diagram with A^ external lines can be generated by a gN<P^ term, so 
that its scaling dimension is [gN\- On the other hand, with only the cjf interaction 
term available, the divergent part of the diagram scales as gYE^. Therefore, we have 



M = [^d^ + ^, eq. (pD. 



If A^ is the maximum power in the potential, the superficial degree of divergence 
is 5 = [fi'Ar](l — V). For the theory to be power-counting renormalizable, it must be 
[qn] > 0, implying 

N<Pd:=^^. (4.7) 

In two dimensions, 6 does not depend on the number of external legs (A^ is uncon- 
strained) and the greater the number of vertices the more convergent is the diagram. 
In four dimensions, the 0^ theory is renormalizable while higher powers of are 
responsible for an infinite number of divergent diagrams. In general, 6 is bounded 
by the dimension of operators which already appear in the bare action. 
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4.2 Actions of fractional field theories 

4.2.1 Role of symmetries 

We are mainly interested in the interplay between the renormalization group flow and 
the symmetries of the model. In ordinary fleld theories, it is natural to impose the 
symmetries of the measure also to the Lagrangian density C These symmetries are 
preserved along the renormalization group flow, and they protect the theory from an 
inflnite multiplicity of non-covariant relevant operators. In fractional fleld theories, 
on the other hand, we are in an unusual situation. The measure does change along the 
RG flow, and fractional Lorentz invariance is not constant with the scale: both the 
fractional charge a and the form of eq. (|2.52| ) flow as well. We do need symmetries 



to protect the model from uncontrollable divergences, but the remarkable role of 
geometry in fractional scenarios strongly constrains the action even in the absence of 
RG-constant symmetries. Calling 11^ and II^ the symmetry group of, respectively, 
the measure and the Lagrangian density, we can construct two inequivalent classes 
of fractional fleld theories: 

• Integer- symmetry scenario. While the symmetry of the measure guarantees 
protection against arbitrary measure distributions, one can prescribe a constant 
symmetry for C in the RG sense. Since in the infrared the Lagrangian should 
be Lorentz invariant, we assume 11/; = Hi, the integer Poincare group in D 
dimensions. 

• Fractional- symmetry scenario. In ordinary models, the symmetry of the mea- 
sure and of the Lagrangian density are the same. If we impose 11^ = 11^, 
we obtain an action invariant under fractional Lorentz transformations. This 
symmetry varies along the RG flow but it forbids non-scalar operators at any 
given a. 

We consider these cases separately but we anticipate that physical arguments in 
favour of the integer-symmetry scenario will be advanced in sections E^ and ff3. 



As in section [ST^, take an action Sa(a) with flxed fractional order at every given 
a. The total Lagrangian can be split into a kinetic and potential part, C = Cjc — Cy. 
By deflnition, the kinetic term is characterized by a dimensionless coupling (in the 
ultraviolet), a quadratic dependence on the fleld 0, and a differential operator /C: 

L^ = U¥:<\>. (4.8) 



Symmetrized kinetic terms of the form D0 D0 can be recast in the form ( f4.8| ) after 
integration by part. Later we shall select the important case 

LK = -\d^<\>d^<\>, (4.9) 
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corresponding to 




(4.10) 



In the last step, the choice oq 7^ a is imphcitly allowed. 
To get a dimensionless operator 



Ok 



fdg^C^, [Ok:]=0. 



the scaling dimension of the field should be 

As far as the potential term is concerned, we take power-law operators 

V J 



where 



[9p] =c?H- 2(^H- [/C]). 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



We will be interested in the left and right fractional Beltrami-Laplace operators 

K,:=r^^''dldZ, IC,:=v"'ood;ood:, (4.15) 

which define an infinite multiplicity of derivative operators: 



^a,7,n 9i 



'^i7i"' 



(ip„ <^(/C^)"0 , 



O, 



Q!,7,n 



9^ 



o,7,n 



(i^a0(/C^)"<^, 



(4.16) 



where a, 7 > 0, 7 will be chosen later, n > 1 is integer, and [5'Q,7,n] = 2(1 — n)'j. 
When n = 1, Oa,'y,i and Oa,'y,i will be denoted as Oq,^ and Oa,'y, respectively. 

4.2.2 Fractional Klein— Gordon equation 

The dynamics of fractional systems have been studied both in Hamiltonian and 
Lagrangian formalism. Fractional phase space and Hamilton equations are known 
for classical-mechanics systems with integer measure and fractional derivatives ||107 - 



[116], with integer derivatives and fractional measure | 117 |, or with measure and 
derivatives both fractional [[118|| - [[12CI|| . Still in fractional mechanics, generalizations 
of the variational principle and Lagrangian equations of motion have been explored 



for one-dimensional actions with integer measure and fractional derivatives |107, 



108| , |114| , |121| , [L22|| , fractional measure and integer derivatives [|117| , |123|| - [[L25|| , and 
measure and derivatives both fractional in one |p.l8| , |119|| and many dimensions ||126|| . 
Variational principle and equations of motion for Lebesgue-Stieltjes actions with 
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absolutely continuous measure were considered in |14, 23, 127 |, in connection with 
fractional dynamics and, in particular, field theory. The same methods allow one to 



construct Noether currents and (non-) conservation laws |T^, |23|, |114| , [I24| , |125|| for 
these dissipative systems [[lO^, |108|, ^, |12|, [128| . 



Variation of fractional actions follow the rules of fractional integration by parts 
(e.g., [ITBl section 2.3.5]). Given two suitable functions /i and /2, from BBl eq. (2.1.50)] 



XO 



XI nxi 

J xn 



(4.17) 



' XQ 



one obtains the integration by parts of fractional integrals of fractional integrands. 
Recalling that left and right derivatives have, respectively, lower terminal Xq = and 
upper terminal Xi = +oo, 

- + 00 



Io,^{f2dV\} 



1 



r(a) 



dxx''-'f2d^h 



1 r+oo 

.«) Jo 



(4.18) 



Using the Leibniz formula, the last derivative term can be expanded as an infinite 
series if 7 7^ 1. In many dimensions and for the scalar theory with potential V{(f)), 
the action integral 

Sa,'y = / d^XVa jC{(j),}C^(f),IC^(j)) (4.19) 



will depend on the kinetic operators /C^ and K.^. Consider an infinitesimal variation 
6(f) which vanishes at the boundary a; = 0, oo. Then, the equation of motion at given 
a 

(4.20) 



_ 0>Jq,7 



6(f) 



IS 



d(J) Va 



dC 



d{}C, 



+ /C, 



dC 



'9(/C, 



(4.21) 



where, for convenience, we expressed the functional variations in terms of differential 
operators of order 27. Varying with respect to operators of order 7 would not change 
the final result but would make its derivation more complicated. 
If 5„,^ = Oa,^ - Oy, eq. (i;^) yields 



= 1 

2 



/C^0 H }C^{Va(, 



r'(0) 



(4.22) 



where V = dV/d(f>. In the particular case 7 = 1, the last two terms of eq. ( [4.21| ) 
collapse into one another and the Euler-Lagrange equation reads 



dC 1 ^ 

(70 Va 



dC 

'°9p0y 



(4.23) 
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where D = dnd^, so that 



= 00 + ^^dy + ^<p - v\<p) 



(4.24a) 



□c 



i^5,0 + (2 - a)(l - «)0 (^ry..^ ) - ^'(0) . (4.24b) 



Here contracted indices are summed over as usual; in particular, the next-to-last 
contribution is ri^^{x'^x'^)~^ = — 1~^ -|- xj"^ -!-■■■ + x]^^. These equations can be 
extended to multi- fractional spacetimes simply by summing over a with weights Qa'- 



= J]^a 






dC 



'9(/C, 



IC^ 



dC 



d{)C, 



(4.25) 



For 7 = 1, eq. ( [4.24|) has a term proportional to □wq, absent in [jT3|, |14|. In the case 

(4.26) 



of ( |4.9| ), in fact, one has 



= /C0 - V'{(f)) = 00 + 



d„v, 



fi'^a 



v'{^) = D^-^—^d,4>-v'{4>) 



x^^ 



and its multi-fractional generalization 




(4.27) 



With respect to [|13], |T^, there is an important difference: the weight Va 



J^ 9^g(x^) is here regarded as a fixed coordinate profile, not a Lorentz scalar. A 
consequence of this fact is that the equations of motion are not Lorentz invariant 
(in a fractional sense, and for any 7) as in the Lebesgue-Stieltjes interpretation of 

Thus, we find ourselves in a 



14 1 , where no particular profile is chosen for Vo 



[X 



situation where the action is invariant under certain symmetries while the equations 
of motion are not. The root of the problem is, of course, the residual f „ dependence 
after integrating by parts. Then, classical physics is not invariant under any of the 
symmetries enjoyed by the action, except in the IR limit where Lorentz invariance 
is restored up to 0(1 — «) terms. ^ 

*One might try to reconcile this situation with symmetry requirements in field theory by regard- 
ing the Lebesgue-Stieltjes formalism of ||l^, Q as a general framework with absolutely continuous 
measures, and fractional models with integer-order kinetic operators as explicit realizations (via a 
particular coordinate presentation of the differential structure) breaking the formal symmetry of 
the general formalism. This would be somewhat analogous to a choice of background metric in 
general relativity: while the mother theory is diffeomorphism invariant, explicit solutions break all 
or most of the symmetries. Similarly, fractional models may be interpreted as explicit realizations of 
a non-standard differential structure. In that case, however, the fractional Lorentz symmetry of the 
fractional measure would be regarded as accidental. Also, fractional models with fractional-order 
kinetic operators, and any of the geometric properties typical of fractional models considered here 
and in p6[ , would hardly stem from a presentation-independent framework. For these reasons, our 
attitude is not to regard fractional theories as a subclass of more general formulations. 
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4.3 Integer-symmetry scenario 



The right derivative is not a derivative operator for the left fractional coordinates 
q^^. Even if one is entitled to define the theory with one type of derivative, the other 
type will always pop in via eq. ( |4.17|) ; this is at variance with the ordinary calculus of 
variations, where the same operator appears both in the action and in the equations 
of motion. 

An alternative is to specialize to the case 7 = 1 in eq. (^]T^) and prescribe that 
the Lagrangian contains only ordinary derivatives, eqs. (f4.23|) -( ^.24b|) . This scenario 
has the double advantage of simplicity (removing any reference to independent sec- 
tors) and of constituting the correct Ansatz for a field theory which is power-counting 
renormalizable at the two-dimensional critical point a*. It is indeed one choice of 
operator d^ against infinitely many others parametrized by 7, but it is also the only 
one with both these characteristics. 

In the spirit of the integer-symmetry class of scenarios, this is also the natural 



choice. Contrary to integer-order actions, however, (f)d^d^(j) 7^ ^9^ 



under in- 



tegration by parts, because the non-trivial measure weight is responsible for extra 
friction terms. Therefore, the ordering of derivatives in the kinetic term is impor- 
tant to determine the dynamical properties of the model (also at the quantum level, 
through the propagator, as we shall see later). Here we take eq. (|4.1CI| ), as in |T3|, |l4| 
but contrary to [|l^ (where the 7 = 1 special case of eq. (|4.15|) , /Ci = D, was 
assumed) . 

The field and coupling ( ^.14| ) have dimension 



2 



[gp] =du-p 



dn 



(4.28) 



The kinetic operator Ojc is marginal. The Op are marginal or relevant if, and only 
if, 

P<Pa:= -r^ ■ (4.29) 



dR-2 



The total action at a given a is 



S, 



a{a) 



dQa{cj) 



^V-'d,^ 



p ^ 



(4.30) 



The fact that the dimension of the scalar field vanishes at rfn = 2, eq. ( [4.28]) , signals 
the presence of a critical point, as it happens, for instance, in Hofava-Lifshitz models. 
For general D, there exists such a point at a = a^, = 2/D, where du = 2 and all 

A critical point, eventually identified with a UV fixed 



the Op are relevant 



H, HI 



point, is important for the power-counting renormalizability of the theory, proven in 
section ^]6] (see also [[l^)- This property supersedes any heuristic argument such as 
eq. (|3.31| ) in selecting dn = 2 as special. 
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In section |3.3| , we have seen that the existence of a norm throughout the dimen- 
sional flow constrains the geometry to configurations with a > 1/2. We can make a 
conjecture and identify the critical point a* with the lowest allowed a; this is most 
natural in a multi-fractional scenario, where one should get continuous access to all 
admissible points in the flow. Then, if the lower limit a = 1/2 coincides with the 
critical value a^, = 2/D, one gets D = 4. Therefore, the dimension of the embedding 
can be fixed to D = A not because of phenomenological reasons, but by virtue of a 
combination of reasonable geometric requirements. 

4.4 Fractional-symmetry scenario 

Instead of fixing the kinetic term to the integer Beltrami-Laplace operator /Ci or to 



one could take fractional differential operators. In doing so, we would accept 
the fact that fractional theories will generically feature different operators at the 
action and the dynamical level, but at the same time we would like to understand 
why. The answer is given by eqs. (|2.7|) - (|2.9|) : The operators d" and 8°' are related to 



each other by a reflection centered at xq + Xi. Debating whether to use d" or 9" in 
the action is tantamount to discriminating, in the integer case, between the operators 
d/dx and —d/dx, and the statement that d"q = f{g) 7^ 1 for some non-trivial /(g) 
is the analogue of dx/d{—x) = — 1 7^ 1. In the fractional case, powers of —1 make 
the issue and the formulae visually complicated, but the essence is the same. In 
this respect, the choice between a left or a right theory is merely political. If the 
integration domain is symmetric, xq = —xi, then also the reflection is symmetric and 
we usually describe it as a parity or a time reversal transformation. Also, the Weyl 
derivative ood"' and the Caputo operator d" with lower terminal Xq = are formally 
conjugate to each other under a reflection at infinity. 



Operators ( 4.16 ) are not invariant under fractional Poincare transformations with 
charge a unless ■y = a. We keep only Ca,a,n, which are based on the kinetic operator 

JCa = v'"'ood;ood:, [£J = 2«. (4.31) 

Then, for du = Da, 

(p=(j-l]a, [(7,,„,„] = -2«(n - 1) , [g^]= fD-p^^^a. (4.32) 

The anisotropic case ao 7^ ca 7^ ctj is not possible for a scalar field, unless one 
introduces dimensionful couplings into the definition of the fractional d'Alembertian. 
These couplings should then appear also in the matrices A(^, acting on the geometric 
coordinates q'^ := (a;^)"f /r(a;^ + 1). We do not consider this multi-scale Lorentz 
symmetry here. 

Equation ( [4.31 ) or ICa are not the only second-order operators invariant under 



fractional Lorentz transformations. In fact, one can also take 



d d 
]C^ = r^^^- (4.33) 
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which has the same scahng properties and yields a very similar model. However, 
in this case one could formulate the theory directly in q coordinates with no need 
of fractional calculus, and the connection with fractal geometry would be somehow 
lost. In the following we do not consider eq. ([4.33|) . 

In the external time picture, the engineering dimension of the scalar field changes 
with cr, thus giving a sort of effective continuous renormalization group flow. In 
general, one should add all possible relevant operators, which would emerge anyway 
at the quantum level. 

From eq. ([4. 321), the only marginal operators are, when a 7^ 0, Oa,a,i = Oa,a 



and Cp^, with p£, given by eq. (|4.7|) , while for a = all the operators are marginal. 
Relevant operators exist only if a 7^ and are the Op with p < p^,. This is the same 
condition as in ordinary D-dimensional field theories: in four dimensions, V ^ (f)"^ 
and V ^ (j)^ are relevant and V ^ (f)'^ is marginal. The total action with fractional 
symmetry at a given a is 



Sa{a) = O 



cx{a) 



P<Pd 



>}C 



a(o-) 



9p 
p 



P<Pd 



(4.34) 



The relevant operators are responsible from making the system flow from the 
UV fixed point. A sharp change in the two-point correlation function of would 
happen when [0] =0; this would typically signal a phase transition across a critical 
point. In turn, a critical point can be naturally identified with the UV fixed point, as 
in Hofava-Lifshitz models. This identification does not guarantee the existence of a 
perturbative UV fixed point (which can be inferred only by explicit loop calculations), 
but it is a positive hint in that direction. 

Apparently, the fractional-symmetry scenario does not have a two-dimensional 
critical point in D > 2. If c?h = Da and D > 2, the phase transition would 
happen only when a = 0, corresponding to Pointland (zero-dimensional manifold, 
no spacetime, no dynamics). This does not imply that there is no non-trivial UV 
fixed point, but it makes its existence less clear. In section ^we shall make a crucial 
extension of the theory such that the limit a — )■ will correspond to a spacetime 
with some residual geometric structure. However, for a < 1/2 there is no natural 
norm and the geometric construction of real- valued a models breaks up progressively 
towards the UV. To summarize, it may be possible to extend the flow down to the 
critical point at a = 0, but one would have to give up to geometric structure anyway. 
This extension may not be even sufficient: Regardless the range of a, the argument 
of power-counting renormalizability fails in the present case, since the maximum 
allowed p = po is finite. Therefore, we have fewer indications about the perturbative 
renormalization properties of the theory. This does not jeopardize the supposed good 
UV behaviour of fractal field models, as we shall see in section W^. 
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4.5 Green function 

The propagator of a real scalar field in Lebesgue-Stieltjes theories with absolutely 
continuous measure was computed in |]13|, |T^ for the kinetic operator ( ^^ but, 



because fractional measures are not Lorentz invariant, not many of the details of 
that calculation fit into the present framework. A cleaner derivation of the Green's 



equation was given in [^ for models with kinetic operator /Ci, but here we shall point 
out some subtleties in relating the Green function with the physical propagator of 
the theory. In order to do so, we briefly repeat the calculation of [|1^ but with the 



kinetic operator /C and the fractional Bessel transform |^ , the correct generalization 
of the ordinary Fourier transform. The final result in momentum space will be the 
same as in |T3|. The steps are the same as in ordinary quantum field theory [ |129[ , 



modulo technical differences, and begin with the partition function. At the end 
we will not obtain the propagator (Green function with causal prescription) but a 
generic Green function for the kinetic operator. This exercise is useful for sketching 
both the momentum structure of the actual propagator (an information sufficient to 
complete the power-counting-renormalizability argument) and the caveats entailed 
in the full derivation of the propagator itself, which will be given elsewhere. 

Consider a real free scalar field with mass m in a fixed-order, isotropic (ao = a) 
fractional spacetime. The action is 

Sa = - dQa{x) (j){x) (/C — w?) (f){x) , (4.35) 

where we omitted the integration domain. The free Lorentzian partition function Zq 
in the presence of a local source J' is 

Zo[J]:= /"[I?0]e^[^"+/'^^"(")^(")'^(")1 =: /" [!?(/)] e'^-^ . (4.36) 



To calculate it, we move to fractional momentum space ^6[. This has the same 
measure ga{k) as configuration space, momenta are non- negative, and there exists 
an invertible transform in terms of Bessel functions of the first kind. Let 



D-l 



Ca{k, x) := r{ao){koxy^^ J^_^{kox^) JJ r(a)(A:V)^-t J|_i(Fx^) (4.37) 



2 



be an eigenfunction of the kinetic operator /C, 

ICca{k,x) = -k^Ca{k,x), k^ := k^kf" = -{kof + \k\'^ . (4.38) 

A reason why to choose /C instead of /Ci, ICa or /C^ is because the momentum-space 
transform is expanded on a basis of c^, which are not eigenfunctions of the other 
kinetic operators. The transform of a function f{x) and the anti-transform are |]36 | 



_ /"+00 r+oo _ 

f[k):= dg^ix) f{x)ca{k,x) , f{x) = dQa{k) f{k)ca{k,x) . (4.39) 

Jo Jo 
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Consistently, the representation of the fractional Dirac distribution in D dimensions 



IS 



C/Qi I i// • i// I 



^Q(y.\'^ } ^ol\'^ 1 ^/Cq( rCj CC j 



(4.40) 



stemming for the definition Sa{x,x') = v~^{x)6{x — x') and the one-dimensional clo- 
sure formula 6{x — x') = x J^ °° dkkJy{kx)Ji,{kx') [ |130| , eq. 1.17.13]. The distribution 
( [4.40| ) acts, indeed, as a delta in fractional space. Transforming both and J in 
eq. ( |4.36| ), one obtains 



S 



J 



I dga{x) / dQa{k) I dga{k') Ca{k,x)ca{k',x) 



1 ~. 



.'2 



(4.40) 



dgaik) 
- / dga{k) 



-^(f){k){k" + m'mk') + J{k)(f){k') 



-^4>ik)ik^ + m^)4>{k) + J{k)~^{k) 



-ipa{k){k^ + m^)ip{k) + 



J{k)J{k) 

f^2 _|_ jy^2 



(4.41) 



where 



^{k) := m - 



^2 _|_ ^ 



(4.42) 



The first term in eq. ( [4.41|) will be a normalization of the partition function. The 
last term can be transformed back to configuration space: 



d0a{k)—r:^ 5- = / dga{k) / dga{x) / dga{x)——- —Ca{k,X)Ca{k,X) 



J^2 _|_ ^ 



J^2 _j_ ^ 

— / dga{x) I dga{x') J{x) Ga{x,x') J{x') , 



(4.43) 



where 



Ga{x,x'; m) :- 



dgaik) 



^2 _|_ ^ 



Cq^ [K^jCj Cq^ [rh^dbj 



(4.44) 



Unlike Green functions in ordinary Minkowski spacetime, Ga does not depend on the 
difference of the coordinates x and x' of the initial and final points. This property. 



unnoticed in [0-[|3 due to the use of a non-factorizable Lebesgue-Stieltjes measure 
and of a non-invertible transform, is a direct inheritance of the measure weight, which 
breaks translation invariance. On a multi-fractional geometry, translation symmetry 
is recovered at large scales [BH]. 
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The free partition function ( 4.36 ) becomes 



ZAJ\ 



\Dlp\ exp 
X exp 



2./0 



r+co 
dQa{x) / dQa{x') J{x)Ga{x,x']m) J{x') 



• /•+00 r+co 

= Zo[0] exp — - / dQa{x) / dQa{x') J{x)Ga{x,x';m) J{x') 
L 2 Jo io 

(4.45) 

The function Ga{x,x';m) obeys the Green equation 

(/C^ — m^) Ga{x, x'; m) = 6a{x, x') , [Ga] = Da — 2 . (4.46) 

The solution given by eq. ([4.44|) is not well defined because we have not specified 
a contour choice. In general, there are infinitely many ways to go around the poles 
and branch cuts of ( f4.44|) in the complex plane (Re/c°, Imfc°), and different contour 
prescriptions correspond to different solutions.^ The most general solution, in fact, is 
a linear combination of two solutions G±. of the homogeneous equation and a particu- 
lar solution G (e.g., the retarded or the advanced propagator) of the inhomogeneous 
equation (f4.4(j|). The arbitrariness of the coefficients of the linear combination corre- 



sponds to the infinitely many possible choices of integration contour. One particular 
choice gives the causal propagator. 

If we had used the kinetic operator /Ci, not only would we have not been able 
to express G^ as a momentum integral, but the left-hand side of the Green equa- 
tion ([4.46|) would not have corresponded to the classical equation of motion ( |4.24|) , 
which has extra terms in the derivatives of the measure weight. Therefore, the usual 
definition of the Green equation as the Klein-Gordon equation in the presence of 
a pointwise source would have no longer been valid. The matching of the Klein- 
Gordon and homogeneous Green equation is related to the issue of microcausality 
(field observable operators commute with one another at spacelike separation, and 
non-local correlations do not give rise to propagation of superluminal messages). In 
ordinary scalar field theory, the combination of the positive- and negative-frequency 
Green functions G+{x — x') = (O|0(x)0(x')|O) and G^{x — x') = (O|0(x')0(x)|O), 
which are solutions to the Klein-Gordon equation, yields the Pauli- Jordan function 
zGpj := G+ — G- = {0\[4>{x) , (j){x')]\0) . The Pauli- Jordan function vanishes outside 
the light cone, thus guaranteeing that the quantum theory is causal. If the func- 
tions G± did not solve the classical equation of motion, the mutual relations among 
different Green functions and their role in causality would have been less transparent. 

Both the contour prescription and causality issues will be reported elsewhere. 
Here, we can nevertheless extract a wealth of physical information from the Green 
function ( [4.44|) . 



An illuminating discussion on the subject can be found in |131 
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Spectrum. The physical spectrum of the theory should be extracted from the 
Feynman propagator, but the pole structure of the Green function, 

G»W--p^. (4.47) 

already points to the final result. Since we are in the half plane Re A;" > 
instead of the full (Re /c°, Im k^) plane, the spectrum has half the usual support 
at k^ = —w?, corresponding to the positive pole 

Re k^ = Vm2+|k|2 . (4.48) 

There is no a-dependence, the momentum-space Green function does not change 
along the dimensional flow, and there is no continuum of massive modes as in 
||14|| .^° Extending to multi- fractional spacetime, the total propagator becomes 
G = Ylia9oiGa or, in the external time picture, G = J dagger) Ga{a)- The spec- 
trum in the integer-simmetry scenario is scale independent. 

Scaling law and critical points. Let A > and consider the scaling transforma- 
tion 

The Green function ( [4.44[ ) transforms as 



Ga (A'^^/'^^x, A'^s/'^^x'; A"'^s/<iH^) = X^^^'^^-^'^'^^G^ix, x'; m) . (4.50) 

This scaling law can be also obtained by the scaling of the fractional diffusion 
equation ( |2.42| ) |TB| and it determines the critical point in the dimensional flow 
at which the Green function is conformally invariant. By definition, it happens 
when 

ds = or du = 2. (4.51) 

In the first case, ao = a = and there is no diffusion at all: this corresponds 
to Pointland, c/h = 0. In the next section we shall extend fractional models in 
a non-trivial way such that even a configuration with a = is endowed with 
a non-singular geometry. Therefore, eventually, one might regard this case as 
physical, albeit its geometry will not possess a norm. 

The second critical point is characterized by a Hausdorff dimension rfn* = 2, 
at the critical value 

a* = I . (4.52) 



^°It is expected that in fractional-symmetry scenarios Ga{k) has an algebraic branch point /c^" = 
— TO^, where the mass coupling has an a-dependence via its scaling dimension. The associated 
branch cut would correspond to a continuum spectrum of massive modes. 
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This is the two-dimensional critical point advertized so far in the integer- 
symmetry scenario, where (is* = 'in* = 2. Taking a^, = 1/2 (the minimum 
allowed value for normed fractional spaces), the topological dimension is con- 
strained to be D = 4. In the fractional-symmetry scenario, this critical point 
either does not exist for D > 2 or it does not correspond to a normed space. 

4.6 Superficial degree of divergence 

Let us make a short summary of the characteristics of multi-fractional quantum field 
theory. First of all, macroscopic physics would not be described by a finite number 
of observables if we included the infinite class of operators Oa,-y,n with fractional 
derivatives, eq. ( |4.16[) . Whether we include these operators or not is a matter of 
definition of the theory, and one is entitled to opt for the formulation with good 
IR behaviour. In the fractional-symmetry scenario there is still a multiplicity of 
operators Oa,a,i, labelled by a, but in the infrared the number of effective couplings 
is finite. In the integer-symmetry scenario there is only one kinetic term throughout 
the flow. Consistently, if one starts from an action with only integer derivatives, 
operators with fractional derivatives never appear. 

Secondly, the renormalization properties of a model are dictated by the dimen- 
sionality of the operator /C in the kinetic term, so that good UV behaviour is guar- 
anteed when the spacetime dimension is the same as the dimensionality of /C. In this 
respect, since the harmonic structure is determined by /C, the spectral and Haus- 
dorff dimensions of the theory (related to, respectively, the harmonic and geometric 
structures ||1^) are equally important and yield complementary informations. 

The scaling argument for the Green function, eqs. ( [4.50|) -( ^.52|) , and the power 



counting of sections [4.3| and ^]4| are in mutual agreement. In fact, the power-counting 
argument of section ^?T] applies, mutatis mutandis, also to fractional theories. One 
difference is in the replacement of the topological dimension D with the Hausdorff 
dimension d^i, due to the non-trivial measure obtained in momentum space: each loop 
integral gives [dga{k)] = Da. The momentum-space propagator has the same scaling 
dimension of the Green function, which is [Go(fc)] = —27 in a general fractional 
scenario with kinetic operator of order 27. In configuration space, G^ contributes 
with a weight Da — 27. For the scalar field theory, interaction vertices do not carry 
dimensionality. Overall, the superficial degree of divergence of a Feynman graph 
with L loops and / internal lines is 

5 = DaL - 27/ < (5max := (Da - 2-f)L . (4.53) 

When a = 1 = 7, one gets the standard result in D dimensions. Otherwise: 

• In the fractional-symmetry scenario (7 = 0), the maximum superficial degree 
of divergence is positive and the power-counting argument is inconclusive re- 
garding the renormalizability of the model. Nevertheless, (5max is smaller than 



-48- 



in ordinary field theory by a factor of a. Using the same hne of reasoning 
leading to eq. ( |4.6D , one can verify that eqs. (^?7|) and ( [4.29| ) are recovered from 

(EH)- 

• In the integer-symmetry scenario (7 = 1), at the critical point a = a* = IjD 
one has 5 < and, at most, logarithmic divergences. If the UV fixed point had 
a < 2/D (non-normed spaces, if D > 4) the theory is super-renormalizable. 

5. Complex fractional theory 

If we regard fractional spacetime models as effective frameworks capturing some 
features of quantum gravity at large, it is important to probe their capabilities in 
that direction, beyond the running of the effective dimension of spacetime. Frac- 
tional models represent continuum spacetimes, and one can conceive applications to 
regimes where the discrete nature of spacetime in quantum gravity models has been 
washed or zoomed away by "hydrodynamical" macroscopic effects. These effects are 
believed to take place when large ensembles of "quanta of space" (spin networks, 
complexes endowed with discrete labels, and so on) are collected together and let 
evolve dynamically. This evolution of a very large number of degrees of freedom 
is presently out of control and scantly explored in most of the theories, but it is a 
promising avenue leading to their yet-unclear continuum limit. Despite their intrin- 
sically continuum structure, can fractional field theories play a role in the effective 
description of this transition? Surprisingly, the answer is Yes. To see this, we need 
a short detour. 

5.1 Prom real to complex fractional order 

A curious feature of the heat kernel trace for a Laplacian on fractals is that it displays 
oscillations. In a metric space of topological dimension D, the return probability at 
small diffusion time a is Weyl's expansion 



V{a) 



(Ana) 



n=l 



(5.1) 



where the coefficients A„ depend on the background metric. From eq. ( p.46|) , do, 



D. For non-smooth sets such as fractals, this expression is drastically modified by 
the presence of discrete symmetries. In particular, for deterministic fractals the 
counterpart of eq. ( ^.1| ) is of the form 

n^) = -^^(^) ' (5-2) 

(47rcr) 2 



where F is a periodic function of Incr |6^, |132|| . Oscillatory behaviour has been 



found analytically and numerically for various fractals [|133|| -[ r3B[| . The phenomenon 
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of logarithmic oscillations | 137 | seems to have two origins. The high symmetry of 
deterministic fractal sets such as diamond fractals and the Sierpinski gasket give rise 
to eigenvalues of the Laplacian with unexpectedly large multiplicity; in turn, these 
are related to the periodicity of the counting measure ||6^, |137| , |138|| . Very recently, 
examples have been found (for instance, the Sierpinski carpet) where log-oscillations 
arise not because of large multiplicities, but because of unexpectedly large gaps in 
the spectrum ||13!J| ]. 

The underlying symmetry mechanism responsible for the oscillations plays a ma- 
jor role in the next development of fractional theory, and we wish to see how it arises 
in that context. We need to recall the relation between fractals and fractional cal- 
culus of real ||140|| - |ll46| and complex order |p.47| , |148|| . This relation was reviewed in 



1^, section 4.4], where the proof of the following theorem was sketched: A fractional 
integral of real order represents either the averaging of a smooth function on a de- 
terministic fractal, or a random fractal support. As a matter of fact, these results 
have been obtained only for fractals embedded in the real line {D = 1) and, to the 
best of our knowledge, there is no literature on multi-dimensional embeddings. We 
do not see any problem in extending the theorem to fractals given by the Cartesian 
product of lower- dimensional fractals, while more general statements might require 
extra work. 

Let J-" be a self-similar set given by A^ similarities (p.9|). It can be expressed 



iteratively as an infinite intersection of pre- fractals: 

oo 

J^ = Pi So-yoSj U) , (5.3) 

k=l k times 

where 

N 

S{U):=[jS,{U) (5.4) 

for any non-empty compact set [/ D J-" such that Si{U) C U. Consider the integral 
of a function f{x) on a self-similar fractal set J-" C [0, 1], 

nx 

Ijr{x) := / dx Vjr[x — x')f{x') . (5.5) 

We temporarily work in dimensionless units ([x] = 0). The kernel vjr depends on the 
geometry of the set and can be determined recursively at any given order of iteration. 
The Laplace transform of eq. ( p.5| ) is 



i^p) := / dxe-P^'lAx) = Mp)fip) ■ (5.6) 

Jo 

Suppose J-" be composed, at the first iteration, by a number of smaller copies of 
length A. The k-th iterate has Laplace-transformed kernel 'O^(p) = Y[n=o9n{p)j for 
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some functions Qn- If all these functions are equal and with argument gnip) = 9{p^^)^ 
and if the asymptotics of g are g{z) ~ 1 + 0{z) for small z and g{z) ~ (71 + 0{z~^) 
for large z (all conditions fulfilled by self-similar and generalized self-similar sets; the 
constant gi is the first probability weight in the self-similar measure ( p.ll|) ), then 
r47|,|l48| 



^fc/-^ .. /-^^ ^-a r. n„ „_N _ l^^l 



lim Vjr{p) = vjr{p) = p °Fa{\np) , a = —— , (5.7) 

fc-i>+oo m A 



where F^ is a log-periodic function [|137|| of period In A: 

.Inp 



Fa{\np + mill X) = Fa{lnp) = y, Qexp j 27rli 

l=-oo ^ 



In A 



(5i 



for some coefficients q. Combining (|5.7| ) with (|5.8| 



27rZ 



Recognizing p~" as the Laplace transform of the fractional weight Va{x) = x"^^/r(a) 
and comparing eq. ( |5.5| ) with (|^), one already sees that Ijr is quite similar to 
a fractional integral J°, were it not for the non-constant contribution ( |5.(j| ). To 
complete the connection, it is sufficient to take the average of the log-periodic function 
Fq, over the period In A: 



/■1/2 
6„ := (F„(lnp)) := / dz F^{ln p + z In X) , (5.10) 



where the value ba depends on the details of g{pX^). The integration range is written 
in a conventional form and it can change according to the choice of units; once the 
log-period is given, the average procedure is uniquely defined. Then, 

{v:f{x)) =ho,—--, (5.11) 

r a 



and 147, 148 



{lTf)= dx'{vAx-x'))f{x') = bJ'^f. (5.12) 

Jo 

This relation is valid for a huge class of sets known as net fractals, and admits 

two other interpretations. Taking the average corresponds to a randomization of 

the fractal structure, where oscillations are cancelled by mutual interference. But 

washing oscillations away can be also seen as dropping all the modes in eqs. (|5.8| ) 

and ( p.9|) except the zero mode / = 0. Now, the cu^ — )■ limit is obtained either as a 

small-similarity-ratio limit, A — )■ 0, or a large-Laplace-momentum limit in eq. ( |5.8|) , 

so that eq. (|5.12| ) can be also regarded as an approximation in Laplace momentum 

space, Ijr ~ 6q,/" as Re(p) -)■ +00 ||140|| - ||146 . 
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Thus, fractional integrals of real order are associated with random fractals. In 
other words, fractional measures either correspond to certain random fractals or, 
alternatively, they approximate Borel measures of self-similar fractals in the limit of 
infinitely refined similarities (continuum approximation), corresponding to neglecting 
the oscillatory structure of fractal kernels. This was the point where the discussion 
in []TB[ ended. However, it is not the end of the story. 

The approximation Ijr ~ (Jjr) ~ /", eq. ( ^.12| ), can be improved by including 
next-to-leading oscillatory modes. Quite generally, integrals on self-similar fractals 
are given by an infinite series of fractional integrals of complex order. Looking at 
eq. ( |5.9|) , the complex fractional measure weight fjr =: Va is naturally defined as a 
sum (or an integral) over frequencies. 



v„[x) 



E 



^u) L/n 



[X 



E 



X 



a— l+io; 



Via + iu) 



(5.13) 



where c^^ are complex coefficients. The zero mode is the real-order measure we 



have considered so far in this paper and in |^, and the average ha is equal to cq. 
Complex measures are obviously not measures in the sense commonly employed by 
physicists: It it not positive (the measure of a set can be a non-negative number) 
and requires a non-trivial extension of the definition of Hausdorff dimension. ^^ To 
explore the properties of these objects, we pick a model with just one pair of conjugate 
frequencies iw: 



,{x) = Co 



X' 



a-1 



r(a] 



+ c. 



X 



a—l+iuj 



„a—l—iuj 



T(a + iu] 



+ 0. 



T{a 



lU) 



(5.14) 



where c^^ = Ic^^le is a complex amplitude, c* is its complex conjugate, and \E' G 



is a phase. We choose to work with eq. (|5.14 ) rather than (|5.13| ) only for simplicity, 
but there may be further justification in doing that. In fact, for Cantor sets the 
three-term weight ( p.l4| ) is a good approximation of the full kernel, where u is the 
average frequency of the leading terms [ 148 1 . 

After some manipulations | p.49| |, one can see that the measure weight ( 15.14] ) is 
real; this happens because one is considering a conjugate combination of complex 
measures. ^^ We set cq = 1 without loss of generality and \E' = 0, commenting on this 
last assumption at the end. Noting that 



Via ± iuj) 



Re 



r(a -|- iuj) 



± ilm 



r(a -|- iu) 



--: Rr{a + iu) ± ilr{a + iu) , 



-'^-'^In this respect, spacetimes associated with such measures are "pre-geometric." We refrain from 
using this adjective because the geometry of complex measures is mathematically well defined, 
although quite different from ordinary geometry. 

^^In general, a self-conjugate real measure has c, 
measure and the ensuing dimension are complex |T3. 



if this condition is not satisfied, the 
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one has 



^n/Ai) \ ^ I 



„a-l 



X 



1 c^e 



iu)\nx 



+ 



CfjC 



-iujlnx 



X 



a-1 



r(a 



X 



a-1 



1 



(5.15) 







r« 



+ 2cujRr(o: + iuj) cosfwlna;) + 2c^lAa + iw) sinfwlna;) 



(5.16) 



The primitive of Va^u) is the oscillatory extension of the measure Qa'- 



where 



A 

^rvAi] . 



X" 



r(a + r 



[1 + Aq,^(^ cos(co'lnx) + -Ba,(^sin(cjlnx)] , (5-17) 



2c,. 



2c.. 



a^ + uj- 



r(a + l)[a-Rr(a + ^w) — C(;/r(Q; + i(jj)] , 
-r(a + l)[u;_Rr(« + ^i^) + a/r(« + i^)] ■ 



(5.18) 
(5.19) 



These expressions are even in w, so we can restrict our attention to positive frequen- 
cies w > 0. In the limit w — ?■ 0, one recovers the power-law measure. Restoring 
dimensionful units, u remains dimensionless, but we must introduce a length scale, 
which we call ^oo? in the arguments of the logarithms. Then, eq. (|5.17|) becomes 



^" 


f x\ f x\ 

1 + Aa^ui COS w In -— + 5a,^ sin w In -— 




=?a,.lXJ r(^^^) 



(5.20) 



5.2 Discrete scale invariance 



Complex fractional models with self-conjugate measure are characterized by oscilla- 
tions governed by a constant 



A. 



exp 



27r 



UJ 



(5.21) 



Notice the highly non-perturbative dependence on the frequency. Asymptotically 

■-!.0+ 



A,,, — -)• i 



-oo , 



(5.22a) 
(5.22b) 



A(^ defines a characteristic (as opposed to fundamental) physical scale as 
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The oscillatory part of the measure weight (|5.16 ) and of the measure ( 5.20 ) is log- 
periodic under the discrete scaling transformation 



. X , X 27m X , , 

In - — — )■ In 1 = In \- nmXu 



u 



n = 0,±l,±2, 



implying 



a; -> A>, n = 0,±l,±2,... 



(5.24) 



(5.25) 



With a slight abuse of terminology, we shall call log-period both the period In A^ and 
the dimensionful scale £^. 

The transformation rule ( |5.25| ) is one of the pivot results of the paper. Log- 
periodicity is a phenomenon intimately related to the presence of a fundamental 
length scale ^oo or? in other words, a microscopic cut-off. This happens due to a 
symmetry unknown to continuous systems or artificially discrete systems such as 
lattices. This symmetry, called discrete scale invariance (DSI), is a dilation transfor- 
mation under integer powers of a preferred, special scaling ratio Aoj ||13?1 , |150| , |151 |. 
Equation ( |5.25|) is a discrete scale symmetry. All deterministic fractals possess a DSI 
by definition p4| , |152|| . For example, the middle-third Cantor set is defined by the 
similarities Si{x) = (l/3)x, S2{x) = (l/3)a; + (2/3). Here, the contraction ratio is 
fixed once and for all: Ai = 1/3. The set is invariant only under contractions with 
ratios A„ = (1/3)", n natural. Log-periodicity and the associated DSI appear also in 
Laplacian growth models, rupture in heterogeneous systems, analysis of earthquakes 
and financial crashes, out-of-equilibrium systems, quenched disordered systems, and 
two-dimensional turbulence ||137|| (for early applications in the theory of phase tran- 
sitions and Levy flights see, respectively, |p.53| , |154|| and ||155|| ). 

Fractional calculus of complex order enjoys (slight modifications of) all the prop- 
erties^^ we have listed in |I6| , section 2.3] but, since the appearance of early papers on 
the subject [|156| , |157| , only recently it has received attention for its physical applica- 
tions ||147|| - |T^ , |158|| . Returning to fractional spacetime models, assume eq. (|5.20| ) as 
the measure along each direction x^. Fractional spacetimes with this measure admit 
two inequivalent interpretations: either as exact per se or as next-to-leading approx- 
imations of a discrete fractal spacetime. In both cases, the extension to fractional 
operators of complex order produces a more complicated geometric pattern unravel- 
ling an underlying discrete scale. One renders the continuous scale invariance ( p.l7| ) 
discrete, and oscillatory fractional measures are much closer to genuine fractals than 
their monotonic counterparts. 

Consistently with ( |5.25| ), to get the continuum limit one should send the fre- 
quency to zero from above, so that the length cut-off vanishes: 

C-> as tj^O+. (5.26) 

"'^'^In particular, left fractional derivatives of order a + iuj with < a < 1 are defined by eq. 
with a replaced hy a + iu! [p8l . 



-54- 



The corresponding asymptotics of the measure weight ( 5.15 ) is 



.0 X 




1 + 2cJ + Lo'c^ 



'^'{a)-ilj^{a) + 2ilj{a)\n 



X 



In 



X 



(5.27) 



where ip{x) = dxT{x)/r{x) is the digamma function. Notice that the measure weight 
at w = is only proportional, and not equal, to its average Va over a log-period. 



5.3 Log-oscillations 

In the case of real fractional order, we have excluded the case a 
measure becomes degenerate fl^, section 2.6], 



Vo[X) 



5{x) 



because the 



(5.2^ 



and spacetime reduces to a structureless point. When the action is equipped with the 
measure ( |5.20|) , however, Pointland is no longer trivial. Setting a = in eq. ( |5.16| ), 
one obtains 



Vo,Lo{Xj 



5ix) + ^ 

X 



Rr{iui) cos ( win -— 1 + Ir{iuj) sin I u\n -— 



(5.29) 



The first term would yield an integration constant (the singular measure of the real-a 
case) and can be ignored, but the rest has a genuine dependence on the coordinate. 

The behaviour of eqs. ( p.l6| ) and ( p. 201 ) is depicted in figure |I[ For u ^ 0, and for 
any a > 0, the weight Va,uj is periodic with increasing period. Due to the power-law 
pre- factor, the amplitude decreases for < a < 1 and is constant for a = 1. On 
the other hand, the measure g^^^ is log-periodic with increasing period and constant 
amplitude for a = and increasing amplitude for a > (inclusive a = 1). The 
amplitudes are magnified by the choice of a large coefficient c^; = 1. In realistic 
fractals, however, c^ is very small and oscillations reduce to tiny ripples around the 
average (see, e.g., ||135|| ). 

It is important to stress a striking difference with respect to real-order multi- 
fractional spacetime. There, like in many other models of high-energy physics, one 
defines a characteristic scale £* (compare with section p.l| ) distinguishing an exotic 
regime ^/£* ^ 1 from a classical regime £/£* <^ 1, depending on the observational 
scale i. Here, on the other hand, these regimes are achieved in a subtler way. The 
characteristic scale i^, determined by the fundamental scale £oo, is the ever-present 
period of the oscillations, and what yields the "classical" result is not an analytic 
expansion in a small quantum parameter, but a spacetime averaging procedure. 

Lifting the assumption \E' = does not present any difficulty. Writing Co := a;+\E', 
one has to replace u with u everywhere (including the characteristic scale icu) except 
in the arguments of -Rr and Jr. The measure is no longer even in u, but one can still 
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Figure 1: The measure weight Va,uj (eq. ( 5.16|) , left-side plots) and the associated measure 
Qa,uj (eq. ( ^-201) , right-side plots) for c^^ = 1, a fixed value of the frequency (here, w = 1), 
and for a = (top row), < a < 1 (center row, with a = 0.95 for Va,uj and a = 0.05 
for ga,u}), and a = 1 (bottom row). The dashed curves are the averages Va = {va,uj) and 
Qa = {Qa,u))- In the case a = 0, the singular term in the measure has been dropped. 



consider positive frequencies. The only non-trivial difference is in the asymptotic 
limit ( |5.27| ), where the trigonometric functions would survive with log-period 27r/\I/. 
Then, it is no longer true that the limit cj — )■ yields a measure proportional to its 
average. This is because the scale £^ does not vanish and the continuum limit is not 
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recovered: eq. ( |5.26|) is replaced by 



^£vt 



>exp 



27r 



as 



w^O" 



(5.30) 



Thus, the interplay of the phase \1' and the zero mode can further enrich the physical 
properties of the model. 

The transient nature of the discrete-to-continuum transition in simplicial quan- 
tum gravity models could be captured by complex fractional models. In fact, the 
lesson from fractal and chaos theory is that, while the large-Laplace- momentum ap- 
proximation and operators of real fractional order are apt to describe "static" fractal 
configurations, log-periodic systems described by complex fractionality also include 
transient phenomena unobservable in the first case. 



5.4 Dimensions 

The parameter a modulates the position and the height of the peaks, but it cannot 
remove them, even in the a = n cases corresponding to Pointland [a = 0) and 
smooth space {a = 1). This reflects the deep relation between genuine discrete 
fractals constructed via contracting maps and their continuum approximation. At 
higher orders in the harmonic expansion, the concept of Hausdorff dimension becomes 
ambiguous. If one used eq. (|2.41|) with the measure ( ^.20|) , the volume V^^\6) of 
a D-hall of radius 6 would strongly depend on 6: even tiny variations of the radius 
would lead to great differences in the output value V^^\5). This value could even 
be negative, and would not monotonically increase with the radius. In general, the 
volume would not be a power law, V^^\6) 7^ 6^^. 

To define the Hausdorff dimension meaningfully, one has to take the average of 
the full measure over a period. Then, dn is simply determined, at all scales, by the 
scaling law of Qa = {Qa,ui)- This is the leading order of the approximation of a highly 
non-trivial fractal measure via an effective fractional measure. In other words, the 
correct operational definition of c?h is not (|2.41|) but 
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ln(V(^)((5)) 



In 6 



(5.31) 



where V^^\S) is calculated with the measure (|5.20| ). Similarly, the spectral dimension 
is the exponent of the leading term in the heat kernel expansion ( p.2| ). In the multi- 
fractional complex case, it is defined through the return probability ( |2.45| ) with the 
measure Qa replaced by the full oscillatory measure in momentum space, and 
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(5.32) 
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This definition replaces eq. ( p.46| ). Again, using eq. (|2.46|) would lead to problematic 



results. The greater the amplitude of the oscillations, the less meaningful the concept 
of spectral dimension would be. One would encounter situations where oscillations 
are so large that neither d^ nor ds make sense any longer (a regime with geometry, 
yet with a very "bad" one), or where ds temporarily becomes greater than da (a 
non- fractal regime) or D. 

To conclude, the Hausdorff and spectral dimensions are the same as in multi- 
fractional spacetimes, because their correct definition entails the average of, respec- 
tively, the volume and the heat kernel over a log-period. 

6. Discussion and research agenda 

In this paper, we completed the construction of a fiat spacetime endowed with cer- 
tain properties typical of fractals. While in [0-[|l5] we outlined the motivations for 



doing so and described some results with general exotic integration measures, in [16 
we focussed on fractional measures and a space with Euclidean signature and fixed 
dimensionality. Here we extended that set-up to Lorentzian signature, changing di- 
mensionality (multi-fractal geometry) and complex fractional order. A scalar field 
theory was given as an example of how fractal geometry deeply modifies the ultra- 
violet structure of a field theory, eventually softening or removing UV divergences. 
For simplicity, gravity was not included in the picture. 

On one hand, fractional field theories can be regarded as effective theories, i.e., 
approximations in the continuum of a microscopic theory either with genuinely fractal 
properties or with a genuine (but non- fractal) dimensional fiow. In this respect, the 
present model would propose itself as a tool to describe effective physics in some 
regimes stemming from fundamental models known to display dimensional flow, such 
as spin foams and others mentioned in the introduction. Then, one would not be 
interested in the renormalization properties of a given fractional action. In this 
case, predictions of fractional models should be associated with features of a given 
full theory, and there would arise the theoretical goal to obtain fractional dynamics 
as an emergent phenomenon. On the other hand, fractional theories may be also 
regarded as fundamental and unrelated to other proposals (fractal or non-fractal), in 
which case one should take care of their UV flniteness. 

This last section is divided in three parts. In the flrst, we describe the hierarchy 
of scales implicit in the measure and the geometric regimes characterizing spacetime 
at different resolutions. In the second and third parts, open issues and applications 
to quantum gravity and non- commutative geometry are discussed. 

6.1 Prom continuum to discrete geometry 

Gathering all the information obtained from the spacetime measure, we can summa- 
rize the emergent physical picture as follows |T^. In complex self-conjugate fractional 
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models, there exists a hierarchy of scales, one fundamental and the others character- 
istic, 

< {L} < L ■ (6.1) 



-oo 



In the simple model with only one frequency, three scales divide six different regimes 
non-perturbatively. We proceed from large to small scales and write expressions in 
one topological dimension, identifying the coordinate x with the resolution i. 

1. Classical regime. At spacetime scales larger than a characteristic scale £*, ordi- 
nary Euclidean/Minkowski geometry and ordinary field theory are recovered. 
The measure along a given direction is 

q{x) ^ {gi^^{x)) = X , i:^i^. (6.2a) 

The number of dimensions can be theoretically constrained to be four. 

2. Multi-fractional regime. At mesoscopic spacetime scales around i^, one obtains 
a model where one can construct field theories on a multi-fractal geometry. The 
scaling of volumes is anomalous and changes with the resolution. The renor- 
malization group properties of field theory on these spacetimes are improved 
in the ultraviolet. The measure is given by 

g{x)^J2^ga,M)), i^L- (6.2b) 

a 

A natural norm for the space exists if the fractional parameter a is comprised 
between 1/2 (at the smallest scales) and 1 (large scales). 

3. Two-dimensional regime. At microscopic scales much larger than a log-period 
£a; but smaller than £*, spacetime is effectively two-dimensional with fractional 
geometry given by the measure 



{^) ~ (^i ..(:^)) ~ x^'^ , L<t:i<L. (6.2c) 



Q 



4. Oscillatory transient regime. In the ultra-microscopic regime i ~ i^^, geo- 
metrical concepts such as dimension and volumes make sense only in average 
(over a log-period), discrete symmetries make their appearance, and despite 
the continuous embedding the notion of continuous spacetime begins to blur. 
In this respect, scenarios with self-conjugate measures are non-perturbative, 
intrinsically quantum models of spacetime. The measure is 

g{x)^J2g^,^{x), ioo<^<L. (6.2d) 
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5. Boundary-effect regime. Here geometry is still given by a fractional continuum 
with discrete scale invariance, but boundary effects become important: this 
happens when the argument x/ioo in the oscillatory part of the measure is of 
order unity. We associate this regime with boundary effects because, according 
to section p.l.3| , a small x expansion corresponds to getting close to the terminal 
at X = 0. Expanding eq. ( p.20[ ) around the point x/ioo = 1, we have 



Qa,ui\X ) 



+ 2Rr{a + iuj) 



1 ^ ^ / ^ 
In -— + O -— 



so that, dropping immaterial constant terms, the measure becomes 

g{x)r^\nx, i ^ ioo ■ (6.2e) 

This is not the same as taking the limit cj — > oo, which is not well defined. 
Notice that all information on the fractional structure of the measure has been 
absorbed in a finite normalization constant. We will comment later on the 
relation of this result with non-commutative geometry. 

6. Fractal regime. Finally, at scales i < ioo, the physics is governed by discrete 
symmetries and the continuum approximation breaks down. From the perspec- 
tive of fractional spacetime, there is no longer a distinction between ambient 
space and its boundary (x ~ 0), meaning that the neighborhood of any point 
X will contain x = 0. This description is clearly inadequate, since the support 
of the measure is highly disconnected (as in dust-type fractals), and the tools 
of fractional calculus must be abandoned definitely: 

g{x)=7, i<ioc. (6.2f) 

Forfeiting some of the above stages, one can even devise scenarios with a shorter 
hierarchy: 



Real-order multi-fractional models. These were the main object of [|16| and 
the first four sections of this paper. The fractional integration order is real 
and one distinguishes between a classical and a multi-fractional regime, via the 
fundamental length £*. In the deep ultraviolet, spacetime is in the continuous 
two-dimensional fractional regime ( |6.2c| ). 

Pure complex models. Setting a = 1, we obtain a complementary scenario 
with just two scales, iuj and £oq. After the oscillatory transient regime, the 
average over a log-period immediately yields the classical result, with no multi- 
fractional structure in between. 
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Both these possibihties have drawbacks. In fractal geometry, one expects both 
anomalous scaling for averaged measures and oscillatory structures. The symme- 
try structure of the real-order multi-fractional model is less rich and it does not 
distinguish different regimes in the UV. Thus, one would loose a number of interest- 
ing connections between some models of quantum spacetime and quantum gravity, as 
we shall comment later. So, the three-scale scenario of the complex multi-fractional 
theory, albeit more complicated, is more complete and, perhaps, better motivated 
than its reductions. 

6.2 Open issues 
6.2.1 Fractal regime 

Fractional operators capture many features of fractals, but not all. In particular, they 
are not complete mathematical realizations of fractals and of diffusion on fractals. 
The features of a genuinely fractal background in the ultra-microscopic regime sym- 
bolized by eq. (|6.2fl) would eventually deviate from those predicted in the continuum 
fractional approximation. For instance, the relation between spectral and Hausdorff 
dimension strongly depends on the class of fractals considered, and it quite often 
happens that dn ^ ds- While in fractal geometry this property is a consequence 
of the definition of the set, in fractional theories it is replicated only in spacetimes 



with anomalous diffusion ||T6[ . In fractional spaces, the type of diffusion depends on 
two choices: that of the operator Pf , entailing a certain degree of arbitrariness, and 
the choice of Laplacian, which is tightly related to the construction of an invertible 
and unitary transform in momentum space. The fractional Bessel transform ( [4.39| ) 
[^ has the desired requisites and selects eq. (|4.10|) as the natural Laplacian. It is 



to be seen whether an invertible transform also exists such that its basis functions 
are eigenf unctions of a fractional Laplacian of order 27. Reflecting upon the inter- 
connected issues of the form of the diffusion equation, the choice of Laplacian, and 
the existence of a well-defined transform will hopefully improve our understanding 
of the transition between fractal and fractional regimes. There do exist fractals for 



which da, = du (e.g., random walks ||159|| and diamond fractals ||135| , |160| ) but, even 
then, real-order fractional models miss the log-oscillations of the heat kernel trace. 
Inclusion of complex modes does improve the approximation and accounts for the 
oscillations. Further exploring the interrelation between complex fractional models 
and general results in fractal spectral theory will benefit our understanding of Fourier 



analysis in both frameworks [16, 161 



In alternative to what done here, one can try a bottom-up approach and attempt 
a brute-force construction on a genuine fractal. Models similar to [^ Q can be 
useful to probe a deterministic-fractal regime, but the technical challenges involved 
therein have partially stalled progress in that direction. Nevertheless, this task may 
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be now within our capabilities, at least in simple scenarios, thanks to the advances 
in fractal geometry p4| , |162|| .^^ 



6.2.2 Scale hierarchy and the Planck length 

The relation between the scales i*, iuj, ^oo and the Planck scale ipi in quantum gravity 
approaches deserves some attention. £^ is a scale derived from ioo via a frequency 
relation, and it acquires infinite multiplicity when considering measures with an 
infinite number of Laurent modes. Therefore, it can be considered as characteristic 
but not fundamental in any obvious sense. This reduces the comparison of ipi to £* 
and ioo- 

Regarded as the numerical constant ( p.30p , the "Planck scale" remains a remote 



concept, since we have not introduced gravity in the picture. Also regarding ipi as a 
symbol for a fundamental scale, the issue is rather undecided. On one hand, from the 
point of view of quantum gravity at large, dimension (in ~ 2 is directly associated 
with the Planck scale (eq. (|3.31|) ), thus preferring the identification of i^, with ipi. 
On the other hand, it is natural to identify 



^Pi 



(6.3) 



both in the perspective of discrete approaches to quantum gravity (because at £oo 
the continuum picture breaks down, while at £* it does not) and in that of non- 
commutative spacetimes (see section |6.3.2| and [§0]). Equation ( |6.3| ) is of interest 



also for the following reason. 

In fractal spectral theory on self-similar fractals, the period In X^ of the os- 
cillations is not an independent parameter, but it is determined by the geometry 
and the harmonic structure of the set [Q (see also [0, section 5.1] for a sketchy 



introduction). Consider a self-similar fractal with measure ( |3.11|) . Instead of defin- 
ing the fractal via eq. ( |5.3| ), one can adopt a formally identical equation where the 
similarities Si are replaced by injections maps fi{x) = riX -|- . . . , where the coeffi- 
cients Ti are called resistance ratios and they are determined by the transformation 
property, in a subcopy i of the fractal, of the Laplacian K, under the mapping /j, 
]C[fi{x)] = r~^K,{x). Summing over all the copies with the appropriate weight, one 
gets lC{x) = Y.i9trtlC[fi{x)] = Y^ilf^fii^)]^ where 

lt--=^ngi. (6.4) 

The harmonic structure is characterized by the probabilities gi of the self-similar 
measure and the resistance ratios Tj, via the combination ( |6.4| ). For the simplest 
self-similar fractals, all 7j are equal. Now, it turns out that the log-period of the 



*In the most general situation, the underlying texture of spacetime may be not just a fractal, but 



a field of fractals tiling the embedding space [163|. For the time being, there is no phenomenological 
reason to consider these configurations. 
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heat kernel trace associated with these sets is, actually, ln7i |Q. So, we can identify 



A^ with 7j, and recognize that the scale hierarchy of fractional spacetimes would 
depend, in a realistic fractal scenario, on the geometric and harmonic structure of 
the underlying set. In particular, one can contemplate the possibihty to obtain the 
Planck length purely from symmetry and geometry. 

6.2.3 Quantum fractional theories 

Dimensional flow was only conjectured in previous papers on fractal spacetimes, due 
to the difficulty in understanding the dimensional properties of these models even at 
a fixed time (here, fixed a). In the introduction, we emphasized the role of dn.s = 2 
in quantum gravity theories; although we have not considered curved fractional sce- 
narios, we see the key principle in action for a scalar, via eq. ( |4.29| ). Thus inspired, 
we built a theory with the aim to obtain a two-dimensional geometry at very small 
scales. In |jT3|, |T^ it was postulated that fractal field theories with ordinary Laplacian 
fiow from a non-trivial UV fixed point at o^h = 2 (corresponding to what we called 
critical point) to an IR fixed point which has d-n = 4 for obvious phenomeno logical 
reasons. Thanks to the requirement of triangle inequality, we have considerably re- 
stricted and better motivated, with respect to |jT3|, |14|, the properties of this fiow. 
On one hand, isotropic fractional geometries with Hausdorff dimension smaller than 
2 can be realized only in embedding spaces with D = 3, which are excluded em- 
pirically. On the other hand, if one assumes that the fiow stops at d^ = 2, then it 
cannot end at dimension greater than D = 4. Embedding spacetimes with D > 5 
cannot reach a two-dimensional UV fixed point. Therefore, if we exclude Kaluza- 
Klein scenarios, the macroscopic dimensionality of spacetime is deeply related to the 
microscopic geometry. Only a full renormalization- group analysis [|164|| -[ p!69|| will 
determine the existence and characteristics of UV fixed points and the perturbative 
renormalizability of multi-fractional models. The results obtained here should allow 
us to begin such a detailed study. 

We have also considered the issue of parity and time reversal in classical fractional 
actions (sections |2.1.1| and [4.4|) . The notion of fractional conjugation [|T6l , 



q^q, a°^9", (6.5) 

is realized simply by an axis reflection from the point of view of the embedding 
coordinates, acting on the space points and its boundary. Just as parity is a dis- 
crete operation rigorously implemented pointwise on a local frame in ordinary curved 
backgrounds, in the general case the conjugation operation will have to be defined 
locally. Since we are not considering curved manifolds, there is a global frame for 
fractional conjugation analogous to the Fermi frame in ordinary Minkowski space- 
time. Conjugation should also play a more involved role at the quantum level, where 
we expect left and right sectors to be mixed in the dynamics. 



-63- 



Other open issues about field theories on fractals should now be prone to direct 
investigation in the easier context of fractional spacetime. For instance, it is not 
known how to place spinors fields on fractals, a problem of interest also in low- 
lacunarity models [^]. The existence of a vielbein formalism would allow one to ask 



the same question even on curved fractional backgrounds. 

Finally, in quantum field theories where Lorentz invariance is broken, pertur- 
bative quantum effects can enhance classical deviations to unacceptable levels. As 
argued on general grounds ||17CI| , |171|| , even if deviations from Lorentz invariance are 



classically negligible, one-loop corrections to the propagator of fields lead to viola- 
tions several orders of magnitude larger than the tree-level estimate, unless the bare 
parameters of the model are fine tuned. For instance, this expectation [p.72|| is indeed 



fulfilled for Lifshitz-type scalar models ||173||. One may ask whether a similar problem 



arises in fractional field theories, but the answer is not obvious to us. Systems with 
a discrete structure may also be protected from large Lorentz violations ||174|| . Even 



if the theory is not integer-Lorentz invariant except in the infrared, the existence of 
other symmetries (discrete first, and then fractional Lorentz) must heavily affect loop 
calculations in a way quite different from a naive modification of dispersion relations. 

6.2.4 Fractional gravity 

The extension of the present framework to fractional curved manifolds is obviously 



desirable for at least the reasons outlined in the introduction and in |jT3|, |T6|, |Tj]. A 
number of properties already displayed by flat fractional spacetimes should be inher- 
ited by fractional models of gravity and tighten potentially interesting relationships 
with theories of quantum gravity. We have seen that a discrete structure of space- 
time emerges through a hierarchy of characteristic scales, and the tools of smooth 
integer geometry become progressively inadequate as the resolution increases. There 
is no obvious obstruction for this to occur also in curved scenarios, barring technical 
difficulties. 



For instance, Tarasov pT] , section 9.2] noticed that failure of the Leibniz rule 



( p.54|) for differentiation in fractional calculus may make the definition of fractional 
differential forms on manifolds problematic. However, this obstacle should be even- 
tually circumvented. In fact, the concepts of parallel transport and Lie derivative 



can be extended to fractional differential calculus [[175|| , and that of manifold exists 



already at the level of pure fractal geometry. Topological spaces such that the neigh- 
borhood of every point is homeomorphic to a neighborhood in a given fractal T are 
called fractafolds ||176| , p.77|| (see also ||178|| ). The existence of fractafolds is a pos- 



itive indication that fractional Minkowski spacetime should admit a generalization 
to manifolds. As far as dynamics is concerned, recasting the Einstein equations in 
fractional fashion [|179|| is not sufficient by itself to define an action theory of grav- 
itation on a fractional manifold, but it should not be difficult to derive them from 
the variational principle. 
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6.2.5 Inflation, big bang, cosmological observations 

A multi-fractional field theory of gravity can have applications in the history of the 
early universe and inflation, which might be explained via the alternative mechanism 
of dimensional flow. In ordinary geometry, the flatness and horizon problems are 
solved if the comoving Hubble horizon contracts while the universe expands at an 
accelerated rate, for a certain period just after the big bang. Inhomogeneities inside 
the horizon are pushed out and diluted by inflation. To realize inflation and explain 
the cosmological perturbation spectra, a slowly rolling scalar field is sufficient. In 
multi-fractional geometry, on the other hand, the effect of a contracting comoving 
horizon can be mimicked by a change in the spacetime geometry, without the need 
to invoke a matter field. The origin of perturbations is less obvious to guess, since 
it would also depend on the gauge symmetries of the model. This scenario will be 
investigated elsewhere. 

Going to smaller scales and earlier times, the log-oscillatory pattern of the mea- 
sure in complex self-conjugate models might shed some light into the big bang prob- 



lem and, hopefully, the Belinsky-Khalatnikov-Lifshitz (BKL) conjecture [ 180 , 181 
Near the classical big bang singularity, a number of cosmological backgrounds admit 
an anisotropic evolution. The simplest BKL model is characterized by three scale 
factors ai(t), i = 1,2,3, which change through a sequence of epochs named after 
Kasner. Within a single Kasner epoch, a^ ~ t^^ and the powers pi obey certain 
conditions dictated by the dynamics. When spatial curvature effects are taken into 
account, one has the following behaviour. Going forward in time from the initial 
singularity, the universe passes through an infinite sequence of Kasner epochs where 
at least one direction is contracting {pi < for one i), although the total spatial vol- 
ume ■\/^ = 0102(33 increases. Across transitions between one epoch and the next, 
the contracting direction exchanges roles with one of the other two, while the third 
evolves monotonically until after a certain number of epochs, constituting a "Kasner 
era." The oscillation amplitudes and the duration of epochs increase during one era. 
The lengths of the Kasner eras are distributed according to stochastic laws which 
can be studied with the methods of chaos theory. In multi-fractional self-conjugate 
theory, we have an oscillatory regime coming neither from the metric tensor nor from 
the dynamics but, still, from geometry. The qualitative behaviour of the measure 
oscillations (increasing amplitude and period) resembles the one described above, 
and the interference of measure oscillations along different directions via an inhomo- 
geneous metric might give rise to something like a BKL pattern near the putative 
singularity. Due to the physically different nature of the oscillatory mechanism and 
to the incompleteness of fractional theory in its present form, it is still premature to 
advance this parallelism any further. 

Finally, observations of the galaxy distribution are generally regarded as com- 
patible with the cosmological principle: at large scales (i.e., i > 10h~^Mpc) the 
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universe is approximately homogeneous and isotropic. However, certain statistical 
analyses of the luminous matter distribution show anomalous correlations at scales 
^ ~ 10 ^ 150/i-^ Mpc |182l-|l88l and, at scales up to 20/1"^ Mpc, matter distribution 



was claimed to have fractal dimension ^ 2 [ 182 |. (Even if anomalous scaling was 
present, by itself it would not mean that matter distribution is the approximation 
of a mathematical fractal. However, in practice, anomalous scaling and fractality 
are taken as synonyms.) At the time of writing, there is no consensus on these re- 
sults, and there is ongoing debate about whether a fractal distribution clashes with 



the cosmological principle or not ||187|| - |1T90| (the upper scale at which homogeneity 



holds has been estimated to be as large as 260/i ^ Mpc ||190|| ). Nevertheless, in case 



of confirmation it would be desirable to have a theoretical model explaining this de- 
viation from the standard lore.^^ The most natural explanation could be found in the 
details of the extreme non-linear regime of gravitational attraction at galaxy-cluster 
scales. However, a cosmological model of fractional gravity should predict deviations 
from standard geometry around scales corresponding to the Hubble horizon at and 
before matter-radiation equality (roughly corresponding to the scales where anoma- 
lous correlations have been allegedly detected), and it could be compared with the 
available data. 

6.3 Other applications 

6.3.1 Doubly special relativity 

Fractional theories are not the first models sporting a non-linear modification of 
Lorentz transformations. If the Planck length is a fundamental building block of 
a theory of quantum gravitation, one may wonder what its significance is in the 
context of special relativity: If it is a minimal length smearing spacetime, should not 
different inertial observers measure the same value £pi? To do so, they should agree 
on an invariant energy/length scale, but ordinary Lorentz transformations act on any 
length-type quantity. With quantum gravity in mind, one can then try to modify 
Lorentz transformations so that the Planck length be observer independent. It turns 
out that the new transformation rules are non-linear in the spacetime coordinates 
and they are parametrized by two invariants: the speed of light and £pi. Frameworks 
implementing the principle of Planck-length invariance are collectively called doubly 
special relativity (DSR) [ p.91| |-[ p!98|] . Global Lorentz invariance in the usual sense 



is only an accidental symmetry of Nature (as first conceived in [ [I99| , |2UU| ) in the 
classical limit £pi — )■ 0. 

If we compare DSR with fractional theories, we notice some common features. 
Lorentz symmetries are deformed and non-linear, and a fundamental length is present 
therein. We have not specified the parameter dependence of the fractional Lorentz 



^^Motivated by earlier large-scale structure analyses invoking a fractal matter distribution, a 
mini-superspace model with varying dimension was proposed in p4|. 
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matrices A^ in eq. (|2.52|) , but they should naturally depend on £00 • In fact, as a 



fundamental building block of the measure, £00 appears in the fractional coordinates 

qi: := QcA^n ■ (6.6) 

At any fixed a and u, the measure can be written in terms of these g^ (rather than 
the coordinates q'^ = (ql^)), and is invariant under eq. ( p.52|) . However, differences 



between doubly special relativity and fractional models soon become apparent. In 
the first case only boosts are deformed, while here all transformations are affected. 
In both cases, Lorentz symmetry is an emergent, non-fundamental property of space- 
time, but in fractional theory symmetries are more structured throughout the scale 
range, from the fundamental DSI of the oscillatory era to fractional Lorentz symme- 
tries of the zero mode of the averaged measure, up to ordinary Lorentz invariance in 
the infrared. Also, while DSR is meant to describe a quantum world, fractional theo- 
ries reach that conclusion from a different route. Fractional geometry automatically 
encodes quantum features, such as a discrete structure at ultra-microscopic scales 
(where DSI emerges) and a natural multi-scale structure equivalently prescribed by 
RG arguments and the lessons of multi-fractal geometry. 

Despite these differences, it may be possible that the DSR paradigm is implicit 
or natural in fractal scenarios. In fact, DSR can arise as a statistical phenomenon 
from a spacetime with a stochastic/fractal structure |f2Ul|| . Further study of fractional 



special relativity, which we have not completely formulated here, will be necessary 
to clarify this quite promising point. 

6.3.2 K-Minkowski 

In D embedding dimensions with integer time direction (ao = 1), the integration 
measure of the action in the boundary-effect regime ( |6.2e| ) is 

dga,oj '~°° VBE{x)d^x := . (6.7) 

Xi ■ --xd-i 

In order to recover this measure in a real-order fractional action with measure Pa, 
one should send a to zero in the measure weight Va and formally keep the leading 
term in the expansion Va{x'^) ~ a/x^, getting 

dga "~ a^'^VBEix) d^x . (6.8) 

However, in the real-order multi-fractional scenario we know that the case a = has 
a pathological geometric structure (a zero-dimensional object): the correct limit in 
the sense of distributions is Va ~ S{x), and the formal inverse-power limit ( |6.8|) is 
at least doubtful. One could simply absorb the vanishing constant a^~^ into a new 
normalization cq for the action, but the geometric considerations leading to eq. ( p.39| ) 
would still lead to tension. 
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The inclusion of the log-oscillations is crucial not only to obtain a finite result, 
but also to give it a completely different physical interpretation. The integration 
measure weight f be can be regarded either as a log-oscillating weight in the boundary- 
effect regime, or as an averaged measure Va = {va,ui) in the limit a — )■ 0. Only the 
first case is well defined. Equation ( |6.7| ) opens up the possibility to link together 
fractional and non-commutative manifolds, in particular K-Minkowski. Imposition 
of a cyclicity- inducing measure in k- Minkowski yields the condition [ 2(J2| 



di[x'v,{x)] = , dMx) = 0. (6.9) 

If one further imposes rotational symmetry, in D — 1 dimensions one obtains v^ix) = 
|x|^^^ ||202|| . However, this is not motivated by strict physical arguments, so another 



solution is Vi^{x) = vbe{x)- 

At this point, it is natural to conjecture a relation between /t-Minkowski and 
fractional models in the boundary-effect regime, with integer time direction. The 
fundamental scale of K-Minkowski (what non-commutativists would call "the Planck 
length") is then identified with i^. Many details should be considered. First, 
eq. ( |3.24| ) implies that the critical point with lowest integer Hausdorff dimension 
in a D = 4 ambient space with integer time has dn = 3. This may seem to be in 
agreement with the fact that the spectral dimension of K-Minkowski is 3 [^, but 
the result of ^ relies on the non-cyclic-invariant measure f^ = 1. One can also try 
to extend the field of investigation and ask whether a mapping exists between gen- 
eral fractional models with power-law measure Qa and non-commutative spacetimes 
with algebras more general than K-Poincare. It turns out that what we called frac- 
tional or geometric coordinates g^ are nothing but coordinates obeying a canonical 
(Heisenberg) algebra onto which non-linear algebras can be mapped. A consequence 
of this identification is that fractional spacetimes are in one-to-one correspondence 
with a certain class of non-commutative spaces. All this is discussed in a separate 
publication [Q. 

6.3.3 Phases of quantum gravity 

Oscillatory measures can have concrete applications in quantum gravity approaches. 
Typically, a quantum spacetime endowed with a fundamental length £* undergoes 



a series of regimes characterized by different spectral and Hausdorff dimensions p3 , 



203|| . Taking the coarse-graining of spacetime texture as physical (i.e., considering a 
quantum non- commutative manifold), the initial condition for the diffusion equation 
of the heat kernel is no longer a Dirac distribution but a Gaussian of width ~ i^, 



^U4| |. In the absence of gravity, the scale-dependent spectral dimension of this object 



is |20l 



^sW = ^^^7^^- (6-10) 
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At large scales/diffusion times, £:»£*, the spectral dimension coincides with the 
topological dimension of spacetime. At scales near the fundamental scale, £ ~ £^,, 
ds ~ -D/2, and if D = 4 one has a two-dimensional regime. At ultra-microscopic 
scales, (is ~ and spacetime dissolves into a zero- dimensional configuration. This 
"trans-Planckian" fuzzy regime was recognized as difficult to interpret in |p03|| . In the 
light of our results, we can provide a description of this regime and it becomes clear 
why we denoted the fundamental scale of |2U3] with £*. Assuming normal diffusion in 
fractional spacetime, the two-dimensional regime corresponds to the (is ~ 2 special 
point in the multi-fractional/RG flow. In a real-order fractional theory, the flow can 
go beyond this point into a no-norm regime, down to the limit a; — )■ 0. In contrast, 
complex fractional theories allow us to pass the artificial barrier at (is ~ 2 and probe 
spacetime at much smaller scales, until the continuum approximation breaks down. 
Estimates of the spectral dimension are not enough to describe the ultra-microscopic 
geometry of spacetime, because they are based on the average return probability, 
eq. ( [5.321 ). Crucially, by looking only at d^ one misses the oscillations of the spectral 
function we do expect in a truly fractal quantum gravity approach. 

As already stated in the introduction and in the previous papers on the sub- 



ject, starting from [|T3[, the logic motivating Lebesgue-Stieltjes field theories is the 
following, (i) At first, one notices that most theories of quantum gravity display 
universal properties of dimensional fiow. (ii) Then, one can ask whether and how 
these properties are related to the problem of UV finiteness. (iii) To this purpose, 
we constructed a geometry and field theory where dimensional flow is an intrinsic 
(not indirect) property, and where one can check the UV finitess explicitly. In this 
respect, fractional theories are independent from other models of quantum gravity 
and forcing two independent models to fit one another (say, fractional theories versus 
quantum Einstein gravity, or versus Hofava-Lifshitz, or versus CDT) may cause mis- 
leading interpretations of otherwise constructive independent insights. Nevertheless, 
it may be possible and very instructive to attempt to establish such a connection. 
With this attitude in mind, in the rest of the section we shall briefiy mention other 
approaches. 

The fractal properties of quantum gravity theories in D dimensions have been 
explored in several contexts. At first, renormalizability of perturbative gravity at 
and near two topological dimensions drew much interest into D = 2 + e models, 
with the hope to understand the = 4 case better | |205| - pTT| . Even in exactly two 
dimensions, there exists a rich fractal structure: short-scale spacetime fiuctuations 
become self-similar, while matter field mass has anomalous scaling governed by an 
order parameter (the "fractal dimension" of spacetime) ||208| , pi2|| . The fractal ge- 
ometry of two-dimensional quantum gravity is difficult to probe in the continuum 
formulation (Liouville theory ||212|| - ||2T7| ), but further progress was made in the dis- 
cretized setting of Z^ = 2 dynamical triangulations | pi8| |-[ p26[] . These studies showed 
the emergence of a branched polymer phase in Z^ = 2 |^1^ , 220 - [P^ , 226 



as well 
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as in D = A Euclidean simplicial gravity [ 227 |- p3T| (dynamical triangulations [232]) 
and CDT [p33| , |234|| (see also [|^). Branched polymers ||235|| are a special case of 
a class of random fractals including tree graphs and random combs. These objects 
have dn = 2, while their spectral dimension is bounded both from above and from 
below, according to the relation 2d}i/{d}i + I) < ds < dji [p36|| - [p4CI|| (see also p41|| ). 
The model of | p242[ is one of the very few instances with scale-dependent spectral 
dimension; the scale hierarchy constructed there bears some resemblance with the 
multi-fractional flow of section |3.2| . 

Studying how fractional calculus encodes these random fractals would help to 
clarify whether a polymeric phase can correspond to scales £00 < ^ < ^* or lower, 
where the continuum approximation breaks down. From the theoretical observations 
we have collected so far, it seems to us that fractional spacetime models might be 
able, by themselves, to provide interesting geometrical insights. Consider, in par- 
ticular, the crumpled phase (phase B) and the two-dimensional branched-polymeric 
phase (phase A) of the CDT phase diagram in 3 + 1 dimensions ||234 |. In phase B, 
the concepts of dimension, metric and volume seem not to play a major role, and 
topology should become important ||230 |. In CDT, a phase-B universe has no exten- 
sion in any direction, ^^ while phase A is characterized by a connected structure with 
(in ~ 2. A third phase (dubbed C) corresponds to a semi-classical four-dimensional 
universe. These regimes coexist in the phase plane and are not sequential in a history 
of quantum evolution. Yet, the regimes of fractional spacetimes can be related to 
phase-space regions as effective "snapshot" descriptions. A comparison with complex 
fractional models indicates that phases A and B of (causal) dynamical triangulations 
correspond, respectively, to the near-boundary regime (formally similar to the struc- 
tureless limit a — )■ in the averaged measure, and where topology should become 
important) and to the oscillatory regime, stuck at the du = 2 fixed point. It is also 
clear why random combs [p39| , p40|, |242|| cannot be associated with phase B: oscil- 
lations are washed away in random structures. The hybrid character of fractional 
theories could make them promising candidates for data fits in causal dynamical 
triangulations. 

Last, the fractional models considered in this paper and in |jl6[ are characterized 
by a Hausdorff dimension dn smaller than or equal to the topological dimension D of 
embedding spacetime, d^ < D. However, there are other scenarios of quantum grav- 
ity where dn > D. An instance is two-dimensional quantum gravity (Liouville theory 
and dynamical triangulations), where typically d-^ > 2 | |219| - p25| (for instance, for 
pure gravity, zero central charge, d^ = 4 ||22U| , |221[ ). Is it possible to accommodate 
these geometries in multi-fractional theories? Technically, the answer is rather sim- 



^^Numerical simulations in 2 + 1 dimensions are visually compatible with the presence of oscil- 
lations in the spectral dimension at small diffusion scales M, but this is a discrete lattice effect 



regarded as unphysical |243]. For the sake of completeness, we mention that, in CDT, fractional 



derivatives have been used as a technical device to calculate the area-to-area propagator [ 244 
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pie. While fractional models such that ao,a < 1 describe spacetimes where D > dn 
(eq. (|3.26|) ), those with fractional charges greater than 1 are characterized by a Haus- 
dorff dimension greater than the topological dimension (including the D = 2 case). 
The details of these models are also governed by the type of diffusion prescribed for 
the fractional manifold which, in turn, dictates the relation between ds and d^. Just 
to give an example, consider a fractional spacetime in D = 2, eq. ( [^.47p : ds = f^dn, 
dn = aQ + a. To get a model with ds = 2 and t/n = 4, mimicking the fractal proper- 
ties of two-dimensional quantum gravity in vacuum, it is sufficient to set ao + a = 4 
and (3 = 1/2 (fractional anomalous diffusion). Such configuration, if desired, could 
be regarded as asymptotic in a multi-fractional setting. This is not sufficient to es- 
tablish a physical equivalence between these fractional models and a particular (and, 
in this case, toy-model) alternative theory of quantum gravity. An in-depth study 
of the relation between fractional theories as effective models and other theories of 
quantum gravity goes beyond the purpose of this work. However, we point out that 
more precise comparative analyses should be well within our capabilities. 
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